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A consistent molecular theory of the non-linear effects of molecular orientation 
produced in dielectric gases and liquids by electric and magnetic fields is given by the 
authors. The molar constant of dielectric polarization Ри, the Cotton-Mouton constant C M» 


the Kerr constant Ka, and the dielectric saturation constant in an electric field Sb and 


in a magnetic field, SM have been calculated for gases (§ § 3 and 4) and for liquids (8 8 5 
and 6). For liquids with molecules possessing axial symmetry, the correlation factors 
Rp, Roy, Ry and Rg associated with the five constants determining molecular interaction 
within the liquid have been obtained. In calculating these factors, no special assumptions 
concerning the nature of the intermolecular forces were introduced. The factors are given 
in the form of functions of the angle Oy between the axes of symmetry of the p-th and 
that of the q-th molecule. In the case of dipole pairwise coupling of the molecules, formulas 
deduced previously by one of the authors and accounting, among others, for the inversion 
of the dielectric saturation ($ 7), are obtained. Finally, a relation connecting the variation 


of the dielectric permittivity in a magnetic field, A в, and the Kerr or Cotton-Mouton 


constant is derived, and the order of magnitude of 4&7, is evaluated. 


I. Introduction 


The present paper is intended to deal with certain phenomena of molecular 
orientation in dielectric liquids subjected to the effect of external electric and magnetic 
fields. We shall consider the following phenomena: dielectric polarization and its 
saturation in an electric and magnetic field, as well as birefringence in an electric and 
magnetic field (the Kerr effect and Cotton-Mouton effect, respectively). 

The theory of dielectric polarisation, the foundations of which, as applied to 
rarefied media, i.e. gases, vapours and dilute solutions, is due to Debye (1912) 
has been generalized for condensed media, 1. е. gases under high pressures, concen- 
trated solutions and liquids by a number of authors, and fundamental papers have 
been published by Debye (1935) and Fowler (1935), Onsager (1936), Kirkwood (1939) 
and Fróhlich (1949). Debye, as well as Fowler, introduce coupling of the molecule 

(209) 
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under consideration with the surrounding medium of a quasi-crystalline type. Denoting 
the energy of coupling by W cos O, wherein @ is the angle between the direction 
momentarily privileged by quasi-crystalline ordering in the environment of the mole- 
cule under consideration within the liquid and that of the dipole, the part of the polari- 
sation due to the dipoles is reduced in the ratio Rp: 1, the reduction factor R, being 


given, for Debye's model, by the relation 


W Et ; 
with y = —, W denoting the energy of coupling, and L — Langevin's function 


kT 
of argument y. 

One of us (Piekara 1938) has shown that, in the case of nitrobenzene, coupling 
with the momentarily nearest molecule plays a greater part than quasi-crystalline 
coupling with the medium (Debye-Fowler coupling). If coupling with the momentarily 
nearest molecule of the same kind tends to produce pairs of almost anti-parallel 
dipoles, then the following relation is obtained for the reduction factor: 


Rp=1— L(y), : (1.2) 


W З 
wherein у = ITO and W denotes the energy of coupling with the nearest molecule. 


This relation has also been proved to be correct by Anzelm (1943). Kirkwood applies 
a macroscopic method to account for short range interactions, and subsequently the 
method was brought to a high degree of generality by Fröhlich. According to their 
theory, the mean polarizability of a molecule (a^? may be written in the form of the 
sum of the induced polarizability and the mean orientational polarizability : 


* 
(at) = af т» (1.3) 
wherein a^ denotes the mean polarizability of a molecule in the condensed-medium 
system, и — its dipole moment when within the medium, and и* — the moment 
induced in a macroscopic sphere by fixing the orientation of one of its molecules. 
If the short range interactions may be neglected, we have u* = и and the Kirkwood- 
Fróhlich theory yields results in agreement with that of Onsager. In the case of A. Pie- 
kara's model of molecular interaction we have 


ие | 
а) =a +3 ap Re» | (1.4) 
with Rp given by (1.2). It is the aim of the present paper to calculate Rp in general 


terms with the fewest possible assumptions as to the mechanism of molecular interaction 
and the nature of the forces. 


In the; second place, the aim of the present paper consists in calculating 
the dielectric saturation, both in an electric and in a magnetic field 


! 


== 
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It will be remembered that the dielectric saturation in an electric field (simply: 
dielectric saturation) resulting from Debye's theory (1939) had been calculated for 
polar gases by Herweg (1920) and subsequently measured in ethyl ether by the same 
author; his results were corroborated by Kautzsch (1928). The effect is to be observed 
as a fall in the dielectric permittivity on placing the medium within a strong electric 
field, the direction of which is parallel to that of the A. C. measuring field. The inverse 
saturation effect consisting in a rise in the dielectric permittivity was found by A. Pie- 
kara and B. Piekara (1936) in pure nitrobenzene. A theoretical explanation of the 
latter, uncommon phenomenon (Piekara 1937a) has been given on the assumption 
already presented here, namely, that the most effective of molecular interactions is 
the coupling between the molecule and its momentarily nearest neighbour. With 
this assumption the following relation is obtained: 


(a) = at +E Rp +3 (97 Rom +2 OF Ек — OF Ree LESS (5) 


4 
и 
3 ee ET : : 
wherein the term — OF Rs = 15 PT Rg dominates and accounts for saturation. 
Assuming for nitrobenzene, as above, that interaction tends to produce approximately 
- anti-parallel pairs of dipoles, we obtain 


Rs =6 — (1 В — £)- (1.6) 
According to this formula, В; changes its sign, becoming negative for y > 1:33, thus 
accounting for the inverse saturation effect increasing the electric permittivity of the 
liquid for sufficiently strong coupling. Recently, the inverse saturation effect has 
been found in certain liquids other than nitrobenzene Piekara and Chetkowski 
1956, Piekara et al. 1957). 

Both theory and experiment prove that the reduction factors, in some cases, 
may take values larger than unity. Therefore it will be more correct to term them 
“correlation factors“, because their appearance in the expressions giving the molar 
constants accounts for spatial and directional correlations of the molecules. 

Generally, for strongly polar liquids the first two terms within the brackets in’ 
(1.5) are small as compared to — Of Rg. The first of them, Of’ Roy, is a purely 
anisotropic term (cf. $ 3 and 6) containing the correlation factor Roy playing an 
important part in magnetic birefringence, i. e. appearing in the molar Cotton-Mouton 
constant. The other, namely 2 0% Rx, is a mixed anisotropic-dipolar term; it contains 
the correlation factor Ёк appearing in electric birefringence, i. e. in the molar Kerr 
constant. Besides these and the purely dipolar term — @$ Rs, the bracket in (1. 5) 
contains deformation terms accounting for deformation of the molecule in the electric . 
field. The necessity of taking electric deformation into account when considering 
dielectric saturation had been previously proved by one of us (Piekara 1937b, 19376). 
Subsequently, electric deformation was considered by Buckingham in a number of 
papers (1956). In the present paper the authors consider both electric and magnetic $ 
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deformation with simultaneous existence of molecular interaction in the condensed 
system. 

The theory of dielectric saturation has been the subject of a number of papers. 
amongst others by Booth (1951) and Schellman! (1957). The latter author restricts 
himself to the purely dipolar term which is definitely calculated on the assumption 
of a specific interaction model, wherein two molecules being the neighbours of the 
one under consideration do not interact with one another. In the theory presented 
in this paper the general formulas for the correlation factors have been obtained 
without making any special assumptions as to the mechanism of molecular interaction; 
solely concerning the form of the molecules it has been assumed that they may be 
considered to be axially symmetric. 

Dielectric saturation in a magnetic field has not yet been shown to exist in true 
liquids?. Efforts in this direction by A. Piekara and M. Scherer (Piekara and Scherer 
1936, Piekara 1936) resulted only in establishing an upper limit of the effect for 
several liquids. Investigations are being carried out with modern means. The theory 
of the effect for paramagnetic gases was given by Van Vleck (1932). The authors' 
theory of dielectric saturation in a magnetic field covers liquid dielectrics, that is, 
diamagnetic substances. Deformation of the molecules in the magnetic field and 
molecular interaction in the dense phase are taken into account. The theory gives 
insight into the connection existing between the effect of dielectric saturation in 
a magnetic field and birefringence in an electric and magnetic field and thus makes 
it possible to predict the effect numerically?. 

Finally, the two remaining effects to be accounted for by the present theory are: 
birefringence in an electric field (the electro-optical Kerr effect) and in a magnetic 
field (the Cotton-Mouton effect). The theory of these effects, based on ideas originating 
in papers by Cotton, Langevin, Debye and Born (Born 1933), applies to gases. However, 
on passing from the diluted system, е. g. a solution of nitrobenzene in a.non-polar 
solvent, to the condensed one, i. e. to the pure liquid, the molar Kerr constant is found 
tod iminish abruptly, whereas the molar Cotton-Mouton constant rises to almost twice its 
previous value (strictly 1.75) (Piekara 1939). One of us has given an expla- nation of 
this behaviour of the molar constants based on the above mentioned model of molecular 
interaction in nitrobenzene (Piekara 1947, 1950, 1951). Here it is the authors’ aim to 
caleulate the molar constants of both birefrigences, accounting for deformation of 
the molecules in the electric and magnetic field and molecular interaction in the 
condensed phase, without introducing any special model. Interesting papers on both 
effects have been published by Buckingham (1956) and Buckingham and Pople (1956). 


! One of us (A. P.) expresses his indebtedness to Dr. Schellman for sending him a typed copy of 
a paper previous to its publication. 0 

2 In the case of liquid crystals the effect was discovered by Jezewski (1924, 1926, 1929) and Kast 
(1924, 1927). 


* Cf. Piekara and Kielich (1957); quite recently similar results have been obtained by op 
(1957). 
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An account of some of the results of the authors’ theory, not comprising, however. 
molecular deformation, has been published in “Le Journal de Physique et le Radium“ 


(1957). 


2. General assumptions of the theory 


Let us consider a given volume V of the dielectric containing N molecules of the 
same species, the positions and orientations of which are given by the variables 
т = T (v, w), v being the positional variables in Cartesian coordinates x, y, z and œ — 
the orientational variables in Euler’s angles 9, Ф, у. The molecules of our system 
possess a permanent electrical dipole moment, and are electrically, magnetically and 
optically anisotropic. 

If the system is subjected to the effect of an external electric field E, each molecule 
undergoes the action of electric forces, besides that of molecular forces. Among the 
latter we distinguish short range and long range forces. In general, the forces of both 
kinds are of an extremely complex character. 

The electric forces comprise the external field E and the field arising from 
polarized molecules in the system. The field considered to consist of the external 
electric field E and of that arising from remote polarised molecules treated as a conti- 
nuous medium is termed the local electric field. The latter is generally anisotropic. 
However, the anisotropy may be considered to recede before the directional inter- 
action of the neighbouring molecules. The local field will be therefore assumed 
isotropic; on the other hand, the treatment of the more important direct mutual 
molecular interaction is entirely free of simplification. 

The molecular forces are themselves strongly dependent on the distances separat- 
ing the molecules and on their spatial orientation. The significance of these forces for 
all effects of molecular orientation in condensed systems, e. g. in liquids, is of import- 
ance; in perfect gases or in diluted solutions they may be neglected, if the coupling 
of dissolved molecules with the solvent is not considerable. - 

The effect of an external electric field is not restricted to producing polarization 
of the molecules as given by the tensor of electric polarizability (а;;); moreover, it 
gives rise to ulterior non-linear electric deformation defined by deformation tensors 


of higher orders (07%) and (с; 


1. 
of intramolecular bonds causing changes in the polarizability and in the permanent 


ial) Such deformation may consist chiefly in a bending 


moment of the molecule. The latter may undergo similar deformations in an external 
magnetic field. Both the electric and magnetic deformations have been taken into 
account in the present paper, as well as the intermolecular forces already mentioned. 

The mathematical treatment of the theory of the effects of molecular orientation 
will be brought under two headings: in the first place, gases and highly diluted 
solutions, in which molecular interaction may be neglected, are considered; in the 
second — condensed systems, composed of mutually interacting molecules, as so- 
lutions or pure liquids. ғ 
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3. Dielectric polarisation and dielectric saturation of a rarefied diamagnetic medium 
in an electric and magnetic field 


Accordingly, we shall first consider a system of N molecules, wherein the existence 
of molecular forces of the mutual interaction of the molecules may be neglected. The 
system is simultaneously subjected to an external electric field E and an external 
magnetic field H making a given angle 2. 

The potential energy of a molecule of the system under consideration in the 
fields F and H may be expressed by the formula 


1 я т 
u(t, F, Н) = — У) ша Е — > У ajaa — > » aj; B; B; Н? — 


1j ij 
iG aN Dijk Gti a a4 P — 5 Y bij. Оч B; Вк Е Н? — 
ijk ijk 
1 ег 1 ет 
-—21 у сум Oi j ® 0 F* — T 3 Cij kl ti €t; Pr Bı Е? Н? E... (3.1) 
ijkl ijkl 


wherein p; denote the components of the permament moment of the molecule in the 
system of reference x, y, z attached to the molecule (i = 1, 2, З), aj; and ар — 
the components of the tensors of electric and magnetic polarisability, respectively, 
We Cie Ofte Сы the components of the tensors of electric and electro- 
-magnetic deformation, and a;, В; — the cosines of the angles formed by the directions 
of the axes x, y, z with those of the fields E and Н; the indices take the values 
И И 

The projection of the total electric moment of the molecule on the direction of 
the field vector E is given by 


or, with regard to (3.1), by 


о unm aj; оң aj Е + = 2 353 о оа 2 Y ma B; p, H? == 


ijk ijk 


1 её ет 
+ (^ Di Cijkl Qi Oj Ope a, ЕЗ + zx Усы a; о; Pr D; ЕН? 3 Є: (3.2) 


ijkl ijkl 


The statistical mean value of the moment mg of a gaseous medium is gna in classical 


statistical mechanics by the formula 


– 5. zm 
ie EL m*.e dc ы; 
DEH pi шор EH) dw (3.3) 


jes dt 
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wherein dt denotes the volume element in positional and orientational coordinate 


> 


зрасе. 


Substitution of (3.1) and (3.2) in (3.3), expansion and integration yield 


(cf. Appendices I and П): 


1 2 
(ME EH = = № С ir a e 


1 


+ (0t" + Of" + Af" + 247") FH? + 
+ (0t + 20% — OF + AT + A AZ) ЁЗ + 


with 
‘ет __ 3 cos? Q— 1 e т e т 
Ө, = 00 ET Ds (За; aij — Qi аз), 
ij 
3 cos? Q—1 m т 2 
ӨТТӨ RT У (Заў ui uj — aii Hi)» 
ij 


'em 1 E v 
Ду" = 30 2 [(3 cos? О — 1) сі"; + (2 — cos? О) ск], 
ij 
lem 1 We 5j 
np — 30 ET 2 [(3 cos? (2 — 1) [555,33 E (2 — cos? Q) Mi 6,2, 
3] { 
ее 1 е е е е 
Өт = 45 kT Ур (Заз; dij — аңа), 
ij 
= E 45 T vm T? à. (3aj; Hi Aj — ай uj), 
e 1 273 
05 45 ЕЗ 45 ЕЗ ТЗ Е Hj ? 
‘= = 9o >) (iij 4-2 cja) 


-amh l (bij; + 2b ji)» 


Q denoting the angle between the vectors E and H. 


(3.4) 


(3.5) 


We define the statistical mean value of the electric polarisability of the molecule 


in the fields E and Н as follows: 


д А 
; (@ DE, = эр (MEDE, H 


(3.6) 
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In particular, the statistical mean value of the electric polarizability of the molecule, 
in weak measuring fields and without biasing fields E and H, is 


{a> = E nis (3.6a) 
The molar constant of dielectric polarizability is, by definition, 
РМ = E Na <at> (3.7) 
with 
e 


being Avogadro's constant, Ny — the number of molecules per unit volume of the 
medium, M — its molecular weight, and d — the density of the medium. 
With regard to definition (3.6a), eq. (3.4) yields 


5 
e»-$Y az) 


whence, with regard to (3.7), the following expression is obtained for the molar 
constant of dielectric polarizability : 


РМ = Sg Y atr e (3.8) 


The molar constant of dielectric saturation in an electric field is defined as follows: 


И, SEE CC ME 


: E (3.9) 


The molar. constant of dielectric saturation in a magnetic field is defined similarly: 


м 4x «аон — (at 
Sem = "gs Na one SD i (3.10) 
These two definitions, together with eq. (3.4) and definition (3.6), yield 
5 = 4 л Мл (Ө? +207 — OF + AY +4 44), (3.11) 
4 Ө." ‘em ү | 
SM = - Na (ӨТ? „КӨТ E Amm Loop am (3.12) 


If variations arising from electrostriction or magnetostriction are not taken into 


account, the electric permittivity in strong fields Е and H is given by the general 
formula 


9 
e(E, Н) =1 +4r № on EDE (3.13) 
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or, with regard to (3.6), by 


OF 
e (E, Н) =1 An № (gn 7 (3.132) 


In particular, with zero field strength, the permittivity is 


oF 
eo = 1 + 47 №, <a?) p (3.14) 


Hence, we obtain the variation of the permittivity arising from the saturation effect 
resulting from a strong electric field applied: 


d wj OF af P 
x stes | 20 | (3.15) 


On the other hand, A ғ", resulting from a strong magnetic field is 


Га: = € (EO) — в == 3 


и 02 8f en 
up ا‎ Е jo 009 


OF 
f=F 5 No 3 (a + (55); zo с) 


It is worth noting that the experimental value of the variation of the permittivity 


A gs EO Eee 


here, f denotes: 


differs from A =, as it represents a superposition of the latter value and of secondary 
effects resulting from electrostriction and the electrocaloric effect; both these second- 
ary effects are relatively small and have been calculated elsewhere (cf. Piekara, Chet- 
kowski and Kielich (1957)). The same is true of the effect of the magnetic field on the 


permittivity. 


4. Refraction and birefringence in electric and: magnetic fields in diamagnetic- 
non-condensed media 


'The component of the total optical moment of the molecule in the direction of 
the field of the light wave E^, if паси іп the fields E and Н аге taken into 
account, is given by: 


mgo = 210 Yi Y; P +25 jk y y; a ЕЕ + 
+32 f nea FF + Fe iu 7:7; Br B, Е Н? 4... (4.1) 


wherein аў, denote the ет of the tensor of optical polarizability of the 
molecule, 65%, cj, and суы — the components of the tensors of electrooptical 
and magnetooptical d ctetu у; —the cosines of the angles between the field 
direction E^ and the axes &, y, z 
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The optical polarizability of the molecule in the fields E and H is defined as 


follows: 


о 9 то (4.2) 


The mean statistical value of a? is given by: 


_ u(t, F,H) 
fare ET POT 


PEE LIA TED EE (4.3) 
I REDE d 


wherein и (т, F, H) is the potential energy of the molecule in the fields E and И, 
as defined by (3. 1). 


By (4.1), (4.2) and (3.1), eq. (4. 3) yields: 
1 о 'om "от 5 
Ga — 3 ) ap 1 (Ot Ar) HET 
+ (Of + OF + AY’ +2 43°) Е? +... (4.4) 
with | 


„т __ 3 cos? Он — 1 у b. Pie 
01 = 90 ЕТ : (3 а а; — а; аў) , 
2] 


, 1 
Ay" = 30 у [(3 cos? Он — 1) сч + (2 — cos? Ән) GL 3 
ij 


toe З cos? Qg— 1 


0 ` (3 aj, а; — ai aj), 
7 
, 3 cos? Qg — 1 o 
ue NOUO У (3 ai; n; uj — ай Hj), 
j 


toe (1 à oe 
Д1 = 30 у [3 cos? Qg —1) су + (2 — cos? Ор) су], 
ij 


Toe 1 oe 
Ay’ = ET M [3 cos? Qg —1) и; bj; + (2 — cos? ОБ) ш; bijl. (4.5) 
1] . "A 


Here, (2р denotes the angle between the light wave vector E’ and the vector E, and Q,,— 
that between the vectors E^ and H. 
The molar constant of refraction is defined by the relation 


RM= Ns pura (4.6) 
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here, (a> is the stasistical mean optical polarizability of the molecule when the 
fields E and Н are absent (ag o ; thus, by (4.4) we may write 
H=0 


4л 1 о 
RM = NA = Da ай. (4.62) 


The molar Kerr constant is defined as follows: 


KM _ 7 Na (Е, 0— «а Ув, о 


: = 47 
and, similarly, the ete Cotton-Mouton constant: 
4. о ; Eas о 
Ces It ү, Ан «аан, (4.8) 


where aj, denotes the optical polarizability of the molecule for a light wave field 
parallel to the vector E or Н (when 2, = 0, or 2y = 0), and «^, — that for a vector 
E^ perpendicular to E ог Н (О; = 90°, ог Он = 90°). 

With (4.4) definitions (4.7) and (4.8) yield 


KM —2 x N.,(O% + Ө + A** +. 2A%), (4.9) 
CM = 2x М, (0?" + At, (4.10) 
where | 


oe 1 о Q5. HH 
0 ES 45 LT Ў (Заза; — aiii) 
ij 


س 


oe 1 о 6 2 
Өз = 45i T? № (За; ui Mj — aii Mj), 


ij 


oe 1 oe ое 
A = 15 X (Biji; — сиу), 
1] 


ое 1 оя ji 
2 — 45kT » (Зи; Dji; — ш biji) 
ij 
aet TET 
01 = 45 LT 3 (Зара; — ара), 
ij 


on T n. om ! 
Ayo = 45 3 (Bij, ij — Cii jj). (4.11) 
m і 


On the other hand, the molar constants КМ and C^ calculated theoretically from 
eqs. (4.9) and (4.10) should;be computed on the basis of experimental results. For 
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this purpose, the relations defining the usual Kerr constant K and the usual Cotton- 


Mouton constant C are employed: 


п (Е, 0) — ny (E,0) -1 


K= с Be 
с_ 210 H)—n,(,H) 1 (4.12) 
n H? 


The light refraction index is given by the fundamental relation 


0 


— 4.13 
H 9 Е 9 ( ) 


n? = 1 + Ал Ny (DE, 


The Lorentz formula may be assumed to describe the local field within the light 


wave: 
2 
Fo CES E (4.14) 
which results in eq. (4.13) yielding the generalized Lorentz-Lorenz formula: 
n?—l 4n | 
ЕТЕР pio Т. E (4.15) 
where 
(az, = } DJ ag (Өт + Ar") H? + 
i 
zi (Ө + Of n Ae en F2 exu (4.16) 


the magnitudes in the brackets being given by (4.5). 


The definitions of the molar constants (4.7) and (4.8) and relations а. 12) апа 
(4.15) yield the well-known relations 


бл? MÍ[EM 
КИ Kaa (z) 
6n? M 


which enable us to compute KM and СМ from experimental data. 


9. The Kerr and Cotton-Mouton effects in a diamagnetic condensed medium 


In a condensed medium, in addition to the potential energy of the N molecules 
of the system in the field of external forces, the potential energy of molecular inte- 
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raction, Uy, should be taken into account. Thus, the total potential energy of the - 
system of N molecules under consideration is equal to 


where 
N 
U (т, F, H) = PU (т, F, H)O (5.12) 


is the potential energy of the system of N molecules in the fields F and H, 
and u(t, F, Н)® — that of the q-th molecule, defined by a formula analogous to (3.1): 


и (т, F, Н)® = — 2i uo F — $ D af; a (9 «@ p 
аворо н... E (5.3) 
ij 


Uy is the potential energy of mutual interaction of the N molecules of the system, 
accounting for intermolecular forces of any kind, e.g. dipolar, inductive, van der 
Waals forces, etc. 

The optical rolarizability 4° of the system is equal to the sum of the optical 
polarizabilities of the molecules: 


N 
A E 9 (5.3) 


where a ?9 is the optical polarizability of the-q-th molecule defined, accordig to 
(4.1) and (4.2), as follows: 


ao — > аё о (0 у® ve > bee, y ye a(? Fu 
+ D a e MO) yf а® " F- 1204 "i y@ yo po 62 H? +. (5.4) 
ji 


In classical statistical mechanics, the statistical mean value (5.3) for a condensed 
medium is given by the following expression: 
U(t, F.H) + UN 


ME к с E 
dion Hcr eT TID TTI Ex 
ORELL AMET Ku dT. 
| with 
| 
dt = П dt, | (5.5a) 


4-1 


dt, — dv, do, is the configurational element; dv, the volume element of the q-th 
molecule is given in the Cartesian coordinates x, y, z; dw, = sin 9, d$, dg, dyp 


| the orientational element given in the Eulerian angles 9, p, y 
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On expanding (5.5) in the field strengths E and Н and on dropping all terms 
containing odd powers of the direction cosines (in isotropic averaging these terms 
vanish, cf. Appendix II), we have 

N 


Г > l r) Г. 
(A SEH = 3 а Oy y jp» тё, Zi. [аза 3 В? 2 © ZR 
1] а=1 Ы] оу: 
N 
ep (a ауен Cy У арар )]) 5 
r=1 s= 


= У (= 2 <. ФУД > + +% 2 < у 3 yp >) ic 
A 
TAN 


N N 
№ у Oya? ee 


а? аы 
“ыйл =i (4 e L < 599050259 У ча Ы 


ijkl q=1 г=1 
E A о, ко 0 oN aj; lk Mp S کک‎ е OO SO 
нту ECÊ У Yonne) +. 
а=1 r=1 s=1 (5.6) 


(X;, i = 1, 2, 3) being the system of refererice at rest attached to the observer, 
(X, i = 1,2, 3) —the molecular system attached to the p-th molecule and (X7, @ = 
= 1, 2, 3) — that attached to the q-th molecule, we have the relation 


gio = == È wg? aP, (5.7) 


where ai^ and a(? are the cosines of the angles between the direction of the field E 
and the axes of the molecular systems attached to the p-th and q-th molecule respective- 
ly, and eo — the cosines of the angles between the axes of both molecular systems 
(Х@, = 1,2,3) and (X,i = 1, 2, 3). Similar relations hold for the remaining 
direction cosines. 

Since the molecules are all of the same kind and in the same conditions, eq. 
(5.6) may be rewritten (for the principal axes) in the form 


и = Маун | (5.8) 


where 


«Ук, н = DE às; n (Orr A; m Н? + 


J- (OT ^* pg, ttp тагара акл (5.9) 
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with 
"жое 3 cos? Qm — || o m a ~ (pq) (Ра) 
9 = 90 kT у а; ау n (За; ОҢ) == )» , 
ij q=1 
"кот 1 от от 
Ay = 30 2 [(3 cos? Он — 1) су» Е (2 — cos? 2н) CH, il 
ij 
toe _ 3 cos? Ор — 1 (ра) cof? D. 
9: 90 kT >, air ay Qu Саз pn 
tape E (a) 09 — S (e? 
өре = I Y aud LEY emo poe. 
ij а=1 r= 
Aye a 30 > [(3 cos? Qg — 1) сй + (2 — cos? Ор) Cn; (5.10) 
nee 


ч sem B У [G cos? Ов — 1) ps bff + (2— cos? Ов) ш BF] 4 af. 


Hence, as in 84, we calculate the molar Kur and Cotton-Mouton constants for 
condensed media: 


K — 2x Ny (Өү* + O;* + Af +243), 


CM = Эл Na (Өрт + At"), 


Т {рг 


1 Ж z 
agers ae од f __ у», 
0 = FT X dij; S (3j Шү Ж 


шү 
with 


ij а=1 
1 N N ( Me 
*oe "deg (р س‎ = o 
87 = ipn 245 23 z d E 
ij 4=1 r=1 


oe 1 ое е | 
Ay iod » (3 сбу — су) 
оё (pa) 
A; == 45 ETT 3 (Зи; bj ju — ш biji) e ex. >: 
| *om __ 1 (ра) on = ў 
i ту; ут айа) $3 ш US 


1 
AT = x я (3c — с). (5.12) 
ў V 
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For molecules possessing axial symmetry, expressions (5. i and (5.12) take 


the form: 

CM — 2л N,(0t" Rom + AP), | (5.13) 
where the terms: 
— aii) (a33 — aii), | 


ое __ 2 (аз 


2 
2 = BPP (азз — а11) и?, 


y 2 
Ay = 45 e + = (3сї1зз + seas]. 
р = 15 zm (b333 "T м, 


9 Ң А, 
1 = BIT (a33 — ау) (а — aii), 


2 


15" = С [сзз,зз + 2cit11 + (361242 — с11,22) + 2 (3с1здәз—— cit33)], (5 14). 


are accompanied by the following correlation factors: 


R= à i: 00 , | k 3 


Ў Aci aene, ES PE 3 
f | TU. EM о vit (ife 3 hé aes ELE ge ten vt UR ME 


ЕлАна 2 
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6. Dielectric polarization and dielectric saturation of condensed media in the 
electric and magnetic field 


The component of the electric moment of the system of N molecules under 
consideration in the direction of the vector E is the sum of the components of the 
electric moments of the molecules: 


N 
о; (6.1) 
а=1 


here, m$? denotes the component of the electric moment of the q-th molecule of the 
system as defined by the formula similar to RS 


E ( (q) .( 
nto Y ща D +}, d ; a9 of @ р ee 5 y bf, a? а! ONT Si 
ij 


ijk 


i 
s b. bem, а® peo pe H? + 


ijk 
1 ) ( 1 (а) (a) pla) p 
Le p сы о af af? а]? Ез + 7 M4 cija 01 aj) BYP B® ЕН? +... (6.1a) 
ijkl ijkl 


In classical statistical mechanics the statistical mean value Mg for a condensed 
medium is given by 
U (7, F, H) + Un 
di f Mge ET dt 
<Мв?вн = ITT AT. (6.2) 
Г e kT dt 
Substituting herein the expressions (5.1) and (6.1), we obtain as in §5: 
(Mig = (A** + O° F + (Өү” + ox -- Ay + 245") Е Н? + 
+ (01* + 203° — Ө; + Дү + 445) Fe +... (6.3) 
With 


МЕ 2 а È a (9 و‎ 


oer Уни 2 Saha 
ө” = ig ahah (УУ a a a a” 


ijkl q=1 r=1 
N N 
i (r) glr) 
-LF Pay Cy wa), 
q=1 pl 
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ijkl q=1 r=1 s—1 
м м М 
«Ууу Quim). 
q=1 r=1 Sem 
1 М 
E ( 
Д = > сы (Ye (4) pg pgs 
ijkl q=1 


ОЛОК e 
a? ај? Ber BP», 


1 м ON 
anms Y AY 
ijkl "HET 
1 и м м 
gt | Jine | > Y a? af? аб P af — 
ОЕ = Ne ёз 
03 У I 
(a) (9) (г) | (v) 
= Xi 0 аһ 01 ) 
Ne a <>. 
1 N N 
Өз = IPT Ay: dij Uk Ш PAS Ps Уа af? a о? al? — 
ijkl 4=1 r=1 s=1 
y N N 
(а) (а) (т) (5) 
-< Xi ipu. spas. 
4=1 2 2 E 
p N № je E N: N 
Өе = У 3 (4) (r) (s) _ (t) 
Sum Hi Hj Lk U1 | c; Qj aE Oe me" 
tjkl Nen uis эй <= t=1 


N 
hy 1 ee (а) (а) „(4 
а, Cijkl Y d aj о хр Ф] 


ijkl q=1 
| N 
de маш Ье = «@ af? а? of? 
6 ET Hi jk! Xj Xk Xì (6.4) 
: ijkl - а=Е`т=1 


Our aim, in analogy to $5, consists in calculating the above expressions 
as functions of the angles subtended by the systems of reference attached to the 
various molecules of the medium (whereas all the angles œ and f are related to the 
vectors E and H). This aim has beén achieved without introducing any assumptions 
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whatsoever as to the nature of the intermolecular forces and without recourse to any 
special model of molecular interaction; we have, however, assumed that the molecules 
may be considered to possess rotational symmetry (an outline of the calculations is 
to be found in Appendix Ш). Denoting by O,,, the angle between the axes of sym- 
metry of the m-th and n-th molecules, we have: 


(MEDE, н = N (MEDE, н (6.5) 
where 
age pe 
CnEPE,Hu = |а 4 ЗАТ Rp| Е + 
+ (9: Rom + OZ" Re + AY" + 243” Ry) ЕН? + 
+ (Of Rom + 207 Rx — OF Rs + AY +447 Rp) ЁЗ +... 59 
with 
lem 3 2 Q е е т т 
Өү" = (a33 — ал) (ass — ат), 
fem ES COS я 
Е = MAME (a33 = ат) Д2, 
Ay” = a; [а + 2 cos? Q) (зеза + 241) +2 (3 cos? © — 1) (олг + 2089) + 
429 (2 — cos? Q) (ста + 2c11,33)], 
2" = тр +2 cos? 0) Ds + 2 (3 cos? 2 — 1) Ыла + 2 (2 — cos? 9) US], 


01 = 


45 BET (аза — ai)”, 


2 2 
Өз = 15 HT (азз — алл) и?, 


4 


O DE 
Өз ~ 45 k3 TS * - 


1 a ee ee 
т = 90 (Зсзззз + 12с113з + 8с1111), 


Ay (b333 + 20113). (6.7) 


= ORT 


The correlation factors Rp, Roy and Ry are given by eqs. (5.15), and 


E NA Que: 9.» у 3 NEE e у ө, cos O, 22! 


т=1 з= 
(6.8) 


P! 
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As was the case in § 3, we employ (6.6) to calculate the molar constant of polarization 


of the condensed medium: 


4л и? 
*М = — e 6.9 


the molar constant of dielectric saturation in the electric field 
Cw = 4л №, (OY Roy +2 OF Re — OF В; + Af +4 A7 Rp) (6.10) 


and the molar constant of dielectric saturation in the magnetic field 


= 4 em em em em 
S = - Na (OF Rem + Of" Rx + 41" + 247" Re). (6.11) 
Similarly are obtained: the permittivity in a weak electric field 
M cL e aw ФЕ 2, 
co = 1 + 4л № (a + зт Re) (sf) (6.12) 
the variation of the permittivity in a strong electric field 
e ey ee TD MER 
A ea; = 3 M PG Qi 2E l — E (6.13) 
and the variation of the permittivity in a strong magnetic field 
LEN p May ОЁ 8f E 
Dlg M о 2E Е — (5 E (6.14) 
where 2 OF 
EM e и LA , 
/=4ам,(а ы) (25). (6.15) 


7. The calculation of the correlation factors in some special cases 


It was seen in $5 and 6 that all five effects due to molecular orientation in con- 
densed dielectrics are described by four correlation factors, which account for the 
mutual interaction of the molecules, namely: 


N 
q= 


Всм = > У, (3 cos? Ө, — D» ; 
q=1 


N N 
Rk = H OI E (3 cos O,, cos O,, — cos On) > 3 
1 М м м x NAUNA N : 
Е => (s SES 2% O 2d pn == О» х А cos Opg cos 64») 


(7.1) 
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where 


UN 


[х(9)е dr 


(X (0)5 = (7.2) 


UN 
Га 
In order to calculate the set of correlation factors (7.1), it is necessary to make specific 
assumptions concerning the form of the molecules and the nature of the forces with 
which they interact. This complex problem is greatly facilitated if an adequate simpli- 


fied model of the mutual interaction of the molecules of the system under consideration 
is adopted. 


The most simple, trivial case is, of course, that of Uy = 0; then 
Кр = km ke k l (7.3) 


and the formulas obtained in § 5 and 6 become identical with those of § 3 and 4 for 
a non-condensed medium, as, for instance, a gas or dilute solution, if the coupling 
with the solvent is to be neglected. 


Another simple case is that of a system consisting of molecules interacting only 
in pairs, each molecule interacting significantly with one nearest neighbour only. In 
this case № = 2 and the set of correlation factors (7.1) takes the form 


Къ = 1 4- соз ОУ; 
3 aM 
Rem = 1 + 5 (eos Эн» — | : 


1 
Re=1+ - (coto — B + 2 (cos Ө), 


1 
Rs =1—3 («о 01,» — i + (4 + 5 (cos Ө) (cos Ө), (7.4) 
with 
| Us 
"Oe ^T dr dr 
(cos" Oj) = Ua HAC риа (1.5) 


n being an integer. In the above formulas 6,5 is the angle between the axes of symmetry 
- of the two molecules, and U, denotes the potential energy of their mutual interaction. 

The exact expression giving the potential energy of mutual interaction of a pair 
of molecules is of a rather intricate nature, and thus it is possible to compute the 
magnitudes in (7.4) and (7.5) only approximately for small values of the energy of 
pairwise interaction — an assumption considerably diminishing the range of appli- 


cations. j 
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A model given formerly by one of us (A. Piekara 1937 a, 1938, 1939a, 1947, 
1950, 1951) involving dipolar interaction of two neighboring molecules, is free of 
this restriction. In this case the interaction energy is that of two dipoles: 


2 
Un = О, = — = (3 cos Ө, cos O, — cos O5) (7.6) 


with 

cos @ = cos O, cos O, + sin ©, sin ©, cos g, (7.7) 
Ө, and Ө, being the angles between the direction of the axis T joining the centers 
of dipoles 1 and 2 and the axes of both dipoles respectively, ф being the azimuth. 


When the dipoles tend to a parallel array one may put ©, = 0, Ө; = Ө (the 
angle between two dipoles) and one obtains cos 0,, = cos Ө and 


О, = — W cos Ө (7.8) 


и? 
with W = 2 — . On the other hand, in the case of nearly anti-parallel coupling we 
r 


put Ө, = Ө, = 90°, ф = 180? — Ө (180? — Ө being the angle between two dipoles) 
and the interaction energy is given by the same expression (7.8) with a different value 
of W and with cos Og = — cos Ө. 


For both cases we have 


(cos" 8,4,» = (+ 1)" L, (у), (7.9) 


where 


п 


f сов" O e? ° sinOdO 


‘Ly (y) 5 ا‎ (7.10) 


л 


f Si HOO 


0 


is a function previously calculated (viz. Piekara 1939a). 
By (7.9), eqs. (7.4) yield the following correlation factors for pairwise coupling: 


L 
Reyes [i552], (7.11) 


T quum WESS 4 Б Г. 


In the above expressions, the upper signs refer to nearly parallel coupling, whereas 
the lower ones correspond to the case of nearly anti-parallel coupling. 
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8. The relationship between A =," and the Kerr and Cotton-Mouton constants 


The variation of the permittivity resulting from the effect of a magnetic field 
and the usual Kerr constant, for a condensed medium, are given by the following 
formulas: 


DONE d= u OF 8f ш 
A Esat = 3 М S H JE | === = эб. ; (8.1) 
уум =з? VS 
жс = KM EASE Еа А МЕ , 
K* = К e xc (8.2) 


On omitting the deformational terms of higher order, the molar constants 57M and 
КМ, in the case of axial symmetric molecules, assume the form (viz. (6.11) and (5.13)): 


SM = ma pp Na (3 cos? Q — 1) д" | д Ком dm Rk): (8.3) 
К*М = Мл д (ôR 4 R (8.4) 
5 5 РМ см + рт RK 


wherein 6”, дапа 6° denoté the magnetic, electrical and optical anisotropies respectiv- 
ely: 


д" = азз — ат, Ó* = a$ — a, 0° = a33 — air ` (8.5) 
The.above equations yield the relations 


SM. 3cos?Q — 1 дт 


=== ем = .6 
em 3 K фо (8 ) 
and 
| a óm 
Деш = (3 cos? Q — 1) (n? +2)? KS E h(s) Н?, (8.7) 
with : 
| 8t miU ESO 7 
h(e) = | == | (5) 2E (8.72) 
and | 


fle) = 4n №, (« V Wr Re) (sf) (8.7b) 


The magnetic anisotropy д” appearing in (8.7) may be determined from the Cotton- 
Mouton constant, which, if the соса term 3” is omitted, has the form 


An 
GM 45 45 LT № дд" Rom. (8.8) 


gj 
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For non-polar substances, it is more convenient to express the relationship (8.7) 
by means of the Cotton-Mouton constant. From (8.3) with u = 0 and (8.8) we have 
immediately 


2() — 
SM c Ssh Е : E CM (8.9) 


which, together with (8.1) and (4.17) yields 


*m 6n? C* д D 
A sq: = (3 cos? Q — 1) ОШО 5 h' (e) Н? (8.10) 
with 
mix of’ 7! OF 
h (s) E [ — | Оё ieu 2E , (8.102) 
oF 
f'(a) === 4л Nya? |25) 2. (8.10b) 


In relationships (8.7) and (8.10) it is necessary to specialize the type of local 
field acting on the molecule. 
For the Lorentz field we have 


e 4-2 
3 


ри Ша» ан = ад, F= E, (8.11) 
where м»; and aĝ; are the values of the components of the permanent moment and 
electric polarizability of the isolated molecule, respectively. 

For the Onsager field we have 


= Зо, Cu б Og, = Lm E (8.12) 
with 
d 2: +1) (ess +2) eT 
{к P ETT CT SEE ss rt 
(2e + 1) (eo + 2) — 2 (e — 1) (e — DA 1 (ао + a02 + боз) 


In the case of isotropic molecules, the degree of anisotropy of the polarizability 
of the molecule A; = 1, resulting in 


ра ES 2 2e +1 
$; = S = | 3 ) eee . (8.13a) 


Hence, for the Lorentz field we have A (e) = 1, whereas for the Onsager field (on 
neglecting anisotropy) 


(26 + воо)? 
h (=) = Зое е) , (8.14) 
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Thus, in the case of polar substances eq. (8.7) may be written for the Lorentz field 
as follows: 
x дт 


A E = (3 coss QESST) (n? t K mm Н? (8.15) 


and, for the Onsager field: 


16 6n (22 + 
{еш = (8 cos 0 — 1) op EE CE at rum E (8.16) 


+ 2)? 8 (2 =? + 22) до 
In the case of non-polar substances, eq. (8.11) assumes the form 
ж бп? є + 2 * дё 
A Esat = (3 cos? О — 1) (n? + 2)? кү С 3 Н? (8.17) 
for both types of local field. 


The ratio of the electrical and optical anisotropies may be evaluated from the 
relationship 


де == 1 


P = maj . | (8.18) 
The evaluation of 4 £7" for nitrobenzene in a field of H = 40 kOe yields 
Aet ore ea 
for the Lorentz field, and 
Aen = 4 x 10-8 


sat 


for the Onsager field. These values differ to some extent from those previously 
given (5. 10-5 and 3 . 10-5 respectively, see Piekara and Kielich 1957) because 
of different experimental date used. 


Appendix I 


Expansion of the statistical mean value of an arbitrary function 


і _ assical statistical mechanics, the statistical mean value of an arbitrary function 
Ф describing the state of a system of molecules and not depending on their kinetic 
energy is given by the following expression: 


Up + UN 
as foe kT dt 11 
< dues UF + UN : (1. ) 


where Up denotes the potential energy of the system in the field of external forces F, 
Uxy— fee energy of mutual interaction of the N molecules of the system, and 


at = H dr, — the volume element of configurational space. 
q=1 # 
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Expanding the exponential functions in power series, eq. (1.1) may be rewritten 
in the form 


(Фук = "= и 
у (ж | ә 
п=0 
with 
U 
[Xe EE dt 
«ХУв-о = «ХУ = oe (1.3) 
e = AT dt 


From (1.2) we obtain, with the desired approximation 


(yr = (Ф) + t a (Ф Ug) — «ФУ QUg)) + 
n (Е ul КФ U2) — 2 46 Ug» (Ug) + (Фу (2 (Ug! (ОЎ + 


l 
21 
soe E) [КФ US) —3 «8 UB) (Ur) +3 ‹Ф Ur) @ Ue — UD) — 
— (®) (6 (Оку? — 6 (Ов) QU£) + QUE)] +... (L4) 
Appendix II 


Calculation of the geometrical mean values of the direction cosines 


If the direction cosines are given in terms of the Eulerian angles 9, p, y, then 
isotropically averaging we have: 


——— 1 А 1, for i=j, 

a; v4 = 3 ду, with д;; = | 0, for ix H (П.Т) 
————— 1 
yj Ar = те (ду ды + Siz ди + ди д»); (11.2) 


E | 
а; а; Pr Bı = 30 [2(2 — cos? 42)д;; ды + (3 cos? Q — 1) (д ди + Ôu 5;x)]; (11.3) 


а; = в; а; X, = а; D; В, = Sex). (IT. 4) 
Similarly, we have: 
1 t 
ind me ОА 
cos ™ = | 2n + 1 ed $ | (П.5) 


0, form —2n +1. 
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Appendix IIT 


Calculation of the correlation factors in the functions 9„„ (the angles between 
the axes of the m-th and n-th molecules) 


Making use of (6.4), we may write for axial symmetric molecules: 


Ave = N isque + (a33 — a11) ) Совт 05 E 
2 
NE СУ cos $5 cos > У 
а=1 
е т т N 
өт" — y (És 4) (аз — ай) | СУ costas cost > — 
2kT эпо 
N 
— (cos? 05» O cos? 23 , 
а=1 
*ет (a33 — a11) и? (X 
Q4" = N == cos? 05 cos 5 cos d» — 
— (cos? 97» (У » cos 9; cos =>) 3 
q=1 r=1 
Өү“ =N (car E Оз cos? 5 cos? 9 T — (соз? 95 > < cos 9 à»). 
N N 
gue COE. и? (< X à. cos? 05 cos 0; cos a» == 
N 
— (cos? 0$» TOPY cos 0; cos >) : 


T N N N A i 
03 = № (s Oy cos $5 cos à x № cos 0, cos я» — 
а=1 т=1 s= 
N N N 
win e D >. cos 05 cos 0; cos O; cos s») з ` (I1) 


а=1 r=1 


E 


where 95 denotes the angle between the axis of symmetry of the p-th molecule and 
the darterion of the electric field E, and 9” — that formed with the direction of the 
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magnetic field Н; similarly we write 07, ди, etc. From the spherical triangle in Fig. 1 


we have 
cos # = cos 05 cos О + sin 05 sin £2. cos p, (111.2) 


Q denoting the angle between the vector E and H. In order to calculate the isotropic 
mean values of the magnitudes (III.1), it is necessary to obtain relations between the 
cosine of 2% (ће position of the p-th molecule being fixed) and the field E, on the one 
hand, and the cosines of the angles 97, ô; and 9; between the q-th, r-th and s-th 
molecules and the field E, on the other hand. From Fig. 2 we have the following 


relations: 


ВЕРХ 


y 
Fig. 1 Fig. 2 
cos 94 = cos 0р cos Ө, + sin 05 sin Ө, cos Pp,» 
cos 97 = cos 0$ cos Ө, + sin 95 sin Oy, cos ф,,, 
cos 9$ = cos 95 cos O,, +- sin 05 sin Ө„ cos Pps- (111.3) 
and 
cos Өг. = cos Ө, cos Op, + sin Ө, зіп O,, cos Py, » 
cos Ө,, = cos O,, cos Ops + sin Op sin O,, cos ф„, 
cos ©,, = cos Ө,, cos Op, + sin Ops sin Ө, cos фи. (1.4) 
By (II.5) and (III. 2, 3, 4), and on averaging, we have 
cos $5 cos 0; = : 0 2 95 дт 1 
р cos ба = > cos Opg , cos 05 = cos 8» = =, 


cos? $5 cos? 0; = se (1 + 2 cos? Ө), 
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е € е 1 
cos? 9$ cos # cos 9, = is (cos Og, +2 cos Op cos Opr), 


a. SS 1 
е е е е 
cos # cos Ù; cos 0, cos 0; = 18 (cos Ор cos O,, + cos Opr cos Өл, + cos Op; cos O,, ), 


TERES edis 
cos? # cos? 9; = 15 [(3 cos? Q — 1) cos? Op, + (2 — cos? Q)], 
(Ш.5) 


т е е 1 
cos? 0р cos à, cos 0; = 15 [(3cos? О — 1) cos Op cos Op, + (2 — cos? О) cos O,,] 


Accordingly, the expressions (II.1) now take the form 


е е 
2 ауу + азз 


A = № 3 : 


он оаа em M TUE Heus Eq A 


тат т A m 2 
Ө" — еа (3 cos? Q — 1) Re, (111.6) 


е е 
азз — а11)? 


Teen m 2( 
Gin NV ЕТ СОМ? 


*ее 2 (азз E aii) и? 
Өт = N gm Re 


*ee pA 
Өз = Мутлу Кз, 


where the correlation factors Rp, Ком, Ёк and К; are given by (5.15) and (6.8) 


КРАТКОЕ СОДЕРЖАНИЕ 

А. Пекара и С. Келих. Нелинейная теория электрической проницаемости и реф- 

ракции диэлектрических жидкостей в электрических и магнитных полях. i 
Авторами дается общая молекулярная теория нелинейных молекулярно- 
‚ -ориентировочных эффектов вызываемых в газах и диэлектрических жидкостях 
приложением электрических и магнитных полей. Вычисляются молярные по- 
стоянные: диэлектрической поляризации Р, Коттона—Мутона См, Керра Км, 
и диэлектрической насыщенности в электрическом SM и B магвитном поле SM для 
газов ($ Зи 4) и жидкостей ($ 5 и 6). Для жидкостей, молекулы которых обладают 
аксиальной симметрией, авторы вычислили факторы корреляции Кр Ксм,Кк,и К S 
относящиеся к этим константам и определяющие взаимодействие молекул 
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в жидкости. В выводе факторов корреляции не сделано каких бы TO ни было 
специальных положений относительно природы сил действующих между моле- 
кулами. Факторы корреляции выражаются функциями угла Opa между осью 
симметрии р-той и а-той молекул. Для случая взаимодействия приводящего к обра- 
зованию пар диполей получаются формулы доказанные раньше одним из авторов. 
которыми учитывается, между прочим, инверсия диэлектрической насыщенности 
(8 7). Наконец, приводится соотношение между изменением диэлектрической про- 
ницаемости в магнитном поле и постоянной Керра или Коттона—Мутона, и дается 


оценка порядка величины Деу. 
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SEMI-MACROSCOPIC TREATMENT OF THE THEORY 
OF NON-LINEAR PHENOMENA IN DIELECTRIC 
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The present paper brings a consistent semi-macroscopic theory of the non-linear 
effects of molecular orientation produced by strong electric or magnetic fields in condensed 
media. The deviation from the quadratic effect of dielectric saturation in a strong electric 
field has been calculated for polar liquids with molecules exhibiting isotropic polarizability 
(5 3). Furthermore, general expressions are derived describing the effect of a strong 
electric or magnetic field on the index of refraction (§ 4) and on the dielectric permittivity 
(§ 5), and yielding the constants of Kerr and Cotton-Mouton as well as the variation of 
the dielectric permittivity in electric and magnetic fields. The molecular interpretation 
of the results yields formulae derived by the molecular method (A. Piekara and S. Kielich). 


1. Introduction 


J. Herweg (1920), applying the theory of Langevin and Debye (1905, 1912) was 
the first to calculate successfully the effect of an electric field on the dielectric constant 
of a dipolar gas. Subsequently, Van Vleck (1932) investigated the effect of electric 
and magnetic fields on the dielectric constant of paramagnetic gases. A. Piekara 
` (1935 a, b; 1937 b, c), moreover, took into account the electric апа magnetic defor- 
mation of the molecule in considering the effect of electric and magnetic fields on the 
dielectric constant of diamagnetic gases. Electric and magnetic deformation of the 
molecule had been considered by M. Born (1933) in his theory of electro-optical and 
magneto-optical effects in polar gases. Recently the problem has been taken up by 
A. D. Buckingham and J. A. Pople (1955, 1956). | 

In investigating the effect of molecular coupling on the dielectric constant and 
on the optical refractive index in polar liquids two tendencies may be distinguished 
as regards the method applied. The first was originated by R. H. Fowler (1935) and 
P. Debye (1935) and developed in papers by Н. Müller (1936), A. Piekara (1937a, 
` 1939, 1950), Н. Friedrich, (1937) A. Peterlin and H. A. Stuart (1939), А. J. Anselm 


(239) 
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(1943) and J. Frenkel (1946), on the basis of a concrete microscopic molecular inter- 
action model and the local field acting upon the molecules. Such a method has 
recently undergone generalization to cover simultaneously polarization, dielectric 
saturation in electric and magnetic fields, and the Kerr and Cotton-Mouton effects 
in polar liquids (Piekara and Kielich, 1957, 1958). 

The other tendency is derived from J. Kirkwood's method (1939) and the much 
more general method of H. Fróhlich (1949), and is based essentially upon the macros- 
copic model of a sphere (containing NV, mutually interacting molecules) immersed in 
a definite continuous medium and subjected to the effect of an external electric 
- field. The method was subsequently applied and developed in papers by A. J. Anselm 
(1944), F: Booth (1951), F. E. Harris and B. J. Alder (1953 a, b), A. D. Buckingham 
and J. A. Pople (1955, 1956) and J. A. Schellman (1957). Both methods apply classical 
statistical mechanics for computing the required mean values. 

The present paper takes into account simultaneously the effect of electric and 
magnetic, linear and non-linear deformation and. of molecular coupling on the 
dielectric constant and refractive index in polar liquids while basing on the latter 
method. From results obtained in this way the author passes to relations obtained 
according to the former method (Piekara and Kielich, 1957). As was the case in all 
above-mentioned papers, the system under consideration is assumed to be in a state 
in which quantum effects are of no importance, thus making the use of classical 
statistical mechanics legitimate. 


2. Fundamental relations 


The dielectric permittivity e of the (homogeneous and isotropic) medium, subjected 


to the effect of a strong electric field E and a magnetic field H forming a given angle 
Q with the former, is calculated from the fundamental equation 


4л 9 (Msg, н _ 4л 9 <Мѕ`еувн 


ATA pex OE EH 
р T WERE Ре spa. Е à 


(2.1) 
where М; denotes the total electric moment of a macroscopic sphere S of volume У; 


and e is the unit vector in the direction of the field E. The statistical mean value 
Ms. е in classical statistical mechanics is given by the expression 


On 
(Ms: e», H=kT x 


JE; la Z (x, Es, Hs _ (2.2) 
where | 
_ UG, Es, Hs) 
Z(t,Es,Hs)=fe № а, (2.3) 


denotes the configurational partition function, Е — the Boltzmann constant, T — the 
absolute température, Es and Н; — the electric and magnetic fields in the presence 
. of the sphere S, respectively. Ù (т, Es, Hs) denotes the total potential energy of 
the sphere 5 in the configuration т with fields Eg and Hg and dt = duda is the con- 


—- 
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figurational element fixing position (volume element dv — dx dy dz) and orientation 
(orientation element dw = sin # 49 do dy in Eulerian angles 9, p, y). Assuming the 
macroscopic sphere S contains Ny molecules, we have 
Ns 
ел. day: dv, = dx 002,5" do, — sin ù dô dg, dy, . (2.3a) 
The effect of a strong homogeneous electric field E on the optical refractive 
index n is calculated similarly: 


4л 9 <М5)в _ 4л 9 < Ms: CDE, 
Vs JE, ED. 9 Ee 


n? — 1 = (2.4) 
oes Му is the total optical moment of a macroscopic sphere S of volume Го, 

— the unit vector in the direction of the field E? of the light wave. . The mean value 
ae ? is obtained from the following relation: 


ЖЕ; д o 
(Ms: e» s = kT —— In Z(r, ES, Es), (2.5) 
dE; 
and 
U(r, ES, Es) 
Z(t, ES, Es)= [ем dr, (2.6) 


where the guantities retain their previous meanings. Similar relations hold for the 
effect of a homogeneous magnetic field H on the refractive index n: in order to obtain 
them it is only necessary to substitute the field H for E in eqs. (2.4) — (2.6). 


3. On the dielectric constant of an isotropically polarizable dense medium in an 


electric field E 


To begin with, let us assume that the homogeneous electric field E acting on the 
system under consideration produces isotropic polarization. This results in the total 
potential energy of the system in configuration t and in the presence of the field E 
taking the form 

U (т, Es) = U (т, 0) — Ms: e Es 5 As ES," (3.1) 
where U(r, 0) denotes the internal potential energy of the system in the absence 
of the external field E (if the sphere S contains № molecules, then U(r, 0) = Ung 
is the potential energy resulting from interaction of the latter), 45 being the scalar 
electric polarizability of the sphere S. 

Substituting (3.1) in (2.3) we obtain 


ASES oo Д 
Z(t Е5) Je us Es 3 R3 ((Ms · e)», (3.2) 
: iU 
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where 


ws) (z, 0) 


(Ms ey = 7 | (Ms:e)'" вт dr, (3.3) 


and 


Ds 2 d 
J= fe (3.4) 


Relation (3.3) gives the mean statistical value of the n-th power of (М; · е) when 
the field E is absent, that is, when all orientations of Mg with respect to е are equally 
probable; hence, on averaging and accounting for isotropic conditions, we may write 


zo ЕА (3.5) 


0 for т = 2n 4-1 


<(Ms е)"> = 


Hence, the configurational partition function (3.2) takes the form 


EI 1 E 
Z (s Es) = Јер: Dr (2л)! Qn +1) de (MS, SO 


and eq. (2.2) yields 


» : Reve 
(Qn +1)! Qn +3) | 
ALTERI ЕЕ. 


E n 
Ness < M3) 


With (2.1) and (3.6) we may write the dielectric permittivity as the following general 
formula: 


Es 2n4-2 
DAC Y TET (s S) (Org 
e—1- y. 3g Ag Es EA . (3.7) 


br d) om 


In the case of a weak electric field, eq. (3.7) yields the following formula for the 


dielectric constant: 
4л. (Ms? OE 
ERS 5 : 
S Vs (45+ ar) E m | (59) 
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With regard to (3.8), eq. (3.7) yields the following expression for the variation 
of the dielectric constant determining the effect of dielectric saturation produced 
within a polar dielectric by a strong DC electric field: 


A o 4 Ес 
тш = سے‎ 


E 30 АТЗ 
—. 70(M$»* — 634M$» (M Mi 
<M) Kari» а 5-901 ы | 3 CLEAR PA 
with 
д О =1 
d(e) == E h- PAR 1 | d (3.92) 
and 
_ 4a (M35 | [AEs | 
Dis y; (4+ ЗЕТ ) EM ; (3.9b) 


here, ¢’ denotes the dielectric constant of the substance in a strong electric field. 


Assuming Fróhlich's model (1949) we have 


4a. MORS За 
(=) е e PIE EUER 


and eq. (3.8) yields Fróhlich's formula for the static dielectric constant’ for small 
field strengths: 


4л Зе «МУ 


о аара A Ee 3.1 
EE. Vote + et ЗЕТ 10) 
and eq. (3.9) take the form 
7 4л 3? Зе AE «М —5 «uy, Bu 
ва = Vs 283 +62, \@в + ғо 30 #3 ТЗ 
2 6 P 2 2.3 | 

pco омо о ыу о ш | вл) 

и 1512 55 T BERS it ы 


where ey, is the high frequency dielectric constant of the medium. 
We now proceed to the microscopic interpretation of the result obtained, Assum- 
ing the sphere S to contain №; molecules, we have 


4s = Sof M. = Yom. (8.19) 
where a? denotes the electric polarizability of the q-th molecule within the sphere 5, 


and MÊ is its electric moment. As all the molecules are of the same species and as 
none of them (except for;a small number on the surface of the sphere) is in any way 
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distinguishable from the rest, we may, in determining the position of an arbitrary 
molecule, write, by (3.12), 
Ns Ns 
(As) = Ns as; (M$) = Ns ну" o» Ope (È È cos Ө" у, (818) 
r=1 5=1 
here, ug is the value of the permanent electric moment of the molecule within the 
sphere S, and O, is the angle between the electric moments of molecules p and 9. 
Finally, by (3.13), we may write eqs. (3.7), (3.8) and (3.9) in the form: 


v = ы-у oh Ene ie pue o cos Ө (Y 3 cos © № 
Za On Ри F3) \ KT м " 
ES ———^— ——— — —— 


r=1 s=1 


(3.14) 
As Es ‘ 
Yes n qx Y mo.) » | 
апа 
ем E 9 Es 
&£ — 1 = 4л a aF ЗЕТ КР Js). (3:15) 
*e e 10 E 
Агш‹ = (Деш) E =; Be RS +. А (3.16) 
where | 
; A Ns | 
(A Esat)ad = = —12л Fs mee d(e) Е? , (3.17) 
and 


Rp = a cos @ me 
| = 


У S » wee e" 


ea 
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ea eod Fue) 


t=1 и= 
NST NS Nise NS 


ау УУ 3 ө, со, + 


r=1 s=1 1=1 y= 


+9 Qi Уу у) Y 3 cos О cos O,s cos в.» | (3.18) 


s=1 ¢=1 и=1 


The relation (3.17) determines the so-called quadratic dielectric saturation effect 
in a polar substance. The second (and following) terms within the brackets in the 
general formula (3.16) determine the departure from the quadratic dielectric satura- 
tion effect as observed experimentally by Kautzsch (1929). The magnitudes Rp, Rs 
and RÉP account for the interaction of the Ns molecules in the dielectric; we denote 
them as correlation factors. | 

For Onsager's (1936) model we obtain 


[£o 12 2е +1 eee Coo 2e +1 
X 3 DPI Fes Quat) sf 3 De £oo x 


кы rapa х= ESD 


2e4+1 ^" с 3(2e2 + £2) 2E l 


where & is the high-frequency dielectric constant of the substance, i. e. 


а (3.20) 
S 


and a, and u, denote the electric polarizability and the permanent electric moment 
of an isolated molecule. Now eqs. (3.15) — (3.17) take the form: 


, Ns 3e Eoo +2 ue Е 
IR = я 3:21 
£ — Sog == АЛ 7 ЕКЕ | 3 | ЗЕТ Кр, ( ) 
and i 
10 3e s: А рођ до E : а) 
obe (KC e Д NOS 0 „е, 22 
Esat = ‘u= daa | Rs L63 =S | 3 | ы da | Ses 
with | 


4 
Е Суз 
з | eT 


. Ns Зе? Se E о 
параша т = 
(26) = H p. = | 


Vs 2e? + $25 
When the induced polarization is. neglected, £% being taken equal to І, eqs. 
(3.21) and (3.22) reduce to 


Ns 36 Мо’ 
4 ры eT BAT Pr 


(3.24) 
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and 
4 t W n EV 
Де SS 3 (A ésar)ad D 63 pem (2) RY E A Й (3.25) 
where 
durs Ns 3e? АА Mes 
а а T gas (seen) aciem P^ Ёз! 


a formula first derived Ьу Van Vleck (1937). 

In the formulas derived previously, the correlation factors Rp, Rs, R® given in 
general terms by eq. (3.18) remain to be calculated in special cases: 

(i) In extremely dilute systems we may assume zero interaction (Uy = 0); then 


Rp=R, = RY = 1 (3.27) 


and eqs. (3.15) — (3.25) hold for polar gases. In this case eqs. (3.21) and (3.24) yield 
the well-known Onsager formula (1936) for the static dielectric constant of a polar 
substance. 

(ii) Concentration is considerable, yet at a given moment the molecules interact 
pairwise only, each molecule having no more than one neighbour in its nearest 
vicinity. Thus, our system is a set of momentary pairs. In this case eqs. (3.18) assume 
the form: 


Rp = 1 + (cos Ө), 
pers | (cos? 9,,» — | + (4 4- 5 (cos Ө,,)) (cos 015), 
RY =1—9 | (сов? OS 2 +3 (4 +7 (cos Ө) (cos Ө» + 


: 
+ zu (70 (cos 6,53 — 63 (cos Ө (cos? Ө + 9 (cos? 6,5). (3.28) 


Thus, our problem reduces to that of computing integrals of the following type: 


Us 
3 f cos" Өззе *T dr, dt ! 
cos Ba curre ECE 


Ù 
f e M dt, dt, 
where the interaction energy of a pair of molecules is of the form 
2 
U, = 00) — Ep (3 cos Ө, cos Ө, — cos Ө) — 


i (3 cos? Ө, + 3 cos? Ө, +2), — (3.292) 


—— 
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and the cosine of the angle formed by the dipole axes of the two neighbouring mole- 
cules is given by 


cos Oy, = cos O, cos Ө, + sin Ө, sin Ө, cos (py — Фф), (3.29b) 


Ө, and O, being the angles between the direction of the axis T joining the centers 
of dipols 1 and 2 and the axes of both dipols, respectively, y = p} — фу being the 
azimuth. 

Because of the involved nature of Uj, integrals (3.29) can be computed only 
approximately (viz. F. E. Harris and B. J. Alder, 1953b, and also A. D. Buckingham 
and J. A. Pople, 1955, 1956). 

In computing (3.29) we apply A. Piekara's method (1937, 1939, 1950) consisting 
in the assumption of a specialized interaction model resulting in a simplification of 
eqs. (3.29 a, b) and hence facilitating computation of (3.29). According to A. Piekara, 
with regard to its chemical structure a molecule may form either a nearly parallel 
or nearly antiparallel momentary pair with its next neighbour; thus, when the dipoles 
tend towards a parallel array опе may assume Ө, = 0, 0, = Ө (the angle between 
two dipoles), which yields cos 0, = cos © and the dipole-dipole interaction energy is 


U, = W,-— W cos Ө (3.30) 


2 
with W = 2 Ix (the permanent dipole-induced dipole interaction is neglected). 
r 


On the other hand, in the case of nearly antiparallel coupling we put 0, = Ө, = 90°, 
ф = 180° — Ө (180° — Ө being the angle between two dipoles) and the dipole- 
. -dipole interaction energy is given by the same expression (3.30) with a different 


value of W and with cos O, = — cos Ө. 
For both cases we have 


cos" Ө) = (+ 1)" Г, (у), (3.31) 


where 


f cos” Ө e» cos Ө зіп O dO À 
І, (у) = = (3.32) 


j e» соз Ө sin © dO 
0 


L, (y) denote functions introduced by A. Piekara (1939b), which are expressed in terms 
VENE W ; d | 
of the Langevin functión L(y), with y — IT denoting the coupling energy of. the 


molecular dipoles in kT units. In general, integration of (3.32) yields - 


ni 1 k еу ر‎ 1-۸-1 AS | 
I, (у) = = |— z) [eec Um rer | (3.32a) 
тл | LES : 
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whence 
1 
= Що P 
2 L 
! P eM Е TS (3.32b) 
H y D 
AERE Se (1-2) ; 
Jd d^ 
By (3.31) and (3.32b), eq. (3.28) may be rewritten as follows: 
Вр=1 Е Г, 


Ro 2(1 =] (GL +4 L, 


x =1—6(1—34] tis GE cr 


+ [sias (15-2) $ (551). (3.33) 
| 2 y yy 


Thus, we-have determined A. Piekara's correlation factors (1937a, 1939b, 1950): 
Rp is the factor corresponding to dielectric polarisation and К; — that corresponding 
to dielectric saturation. In the above expressions, the upper signs refer to nearly 
parallel coupling, whereas the lower ones correspond to the case of nearly antiparallel 
coupling. A. Piekara has shown that the correlation factors Rp and Rg account well 
for experimental data: this is especially important in the case of the factor К $, which 
changes its sign when the interaction energy of a pair of almost antiparallel dipoles 
assumes a critical value, thus accounting theoretically for the positive dielectric 
saturation effect resulting from experiment. 


4. On the refractive index of a dense medium subjected to deformation in an electric 
or magnetic field 


A material medium undergoing the effect of a light wave E^ and that of a strong 
electric field E forming the angle Q, with the vector E^ is subjected to electro-optical 
deformation (apart from mechanical deformations). The total potential energy of the 
system in the state described above may be expressed (in tensor notation) as follows: 


U (т, E$, Es) = U (t, 0) "rcl (42, ES E? + AEE) + 
л | 
— gj (Bory Ёз Е: Е, +3 Bory Eg Ez Ey +...) + 


1 оо 0 то po Р ое 
D 141 (Comme Eo E; E, Eo +6 Corre EE EEE; == „ө ) + *++ 9 (4.1) 
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where U(t, 0) denotes the internal potential energy of the system, M, — the electric 
moment of the macroscopic sphere S, 4‘, and Ag, — the respective electric and 
optical polarizability tensors of the sphere S, В and Сото — the tensors of electro- 
optical non-lineal deformation of the sphere S; the summation indices take the 
values 1, 2,3. М. A‘, AS, Ве апа Сот are tensors symmetrical in all indices. 

On substituting (4.1) in (2.6), expanding in powers of E$ and Es and isotropically 


averaging, we obtain 
д е u ج‎ 2 о ро? 
Т A E; == 4@› Е; Е 


1 
АТ (2 — cos? 26) der б + (3 cos? Ов — 1) (dor дь 4-04,0,)] x (42) 


Z(t, їй Е,) == 


ое 1 с e oe As. M, M, 
X hee ot kT (4: Ave ae 2Воть M, xis Ur Ey ы Aint Ma) | z СНЕ l.. ` 
where 


mo lforo-—r 4 
^ lO for o Tt E 


By (4.2), eqs. (2.5) and (2.4) yield the general formula for the optical refractive 
index of a dense medium: 


| 4л 9 ES |1 я EXS E bris 
2m = — — —€— Е 
АА 1 Vs o Eo Е «Ав (1 SI LT (An TE М, 2 


mI s e — cos? Ок) Ee E rm ; (ans. M, == Ase Ai F ey + 


+ (3 cos? Qg — 1) te om ЫП Mr + А. Aor + 8 a а +. |. 
(4.4) 


Similarly, a formula accounting for the effect of a magnetic field on the refractive 
index may be obtained; in particular, for a diamagnetic dense medium we have: 


re T ә Eo E «Ав (1 = 3117 47 | E 
uS $ < @—соз*2н) ел. T zl " 


+ (3 cos? Он — 1) | „+ Sere) > a h (4.5) 


where Am -is the diamagnetic polarizability tensor of a diamagnetic sphere S, 


and Сот, is that of the magneto-optical non-lineal deformation of the latter; 
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Qy denotes the angle formed by the light wave vector E^ and the magnetic field 
vector H. The tensor C% „ will be symmetrical in the separate pairs of indices 
o, t and m © 


With the sphere S containing Ng molecules of the same species, we may write 


N Ns 
M,— Yue, a У а, 
g=1 а=1 
Ns Ns 
DE (46) 
q= а= 


here, uw, a®, 6%, and с@, are the respective tensors referring to the q-th molecule 
within the macroscopic sphere S. 


Let us attach a molecular coordinate system (X;, 1 = 1, 2, 3) to each molecule 
within the sphere S, and let us transform the tensor components (4. 6) from the external 


(X, о = 1, 2, 3) to the molecular system according to the transformation formula 


Толо = 2) Og; 0,5 O yy Op (m (4.7) 
1] ‘ 

with the assumption that the systems (X,) and (X;) are rectangular, o, c»; 

representing the cosines of the angles formed by the axes X, and X;. By (4. 6) аң 

(4.7) and proceeding as in the case of eq. (3.13), we may rewrite (4.4) and (4.5) as 

follows (assuming the axes of the molecular system attached to the molecule coincide 


with the principle axes of the latter): 


n*—1]1 - 4л NsJ1 : ۴ У off 
ED DI n 2; ай + 30 оу, [в bp nu (as is Qu 


Ns Ns 


2 у 2а 
qi аќ x P4 <> Ba? ogo E ES и vA OF Д (30020 о? — "m. а 


+ (2 — cos? Ор) (as ч 2 2) o>) Esa ar (4.8) 


and 


n2— 1 4л Ns} 1 о 1 Qus 
жүз Раз às. [eem 
1] 


1 


Ns. 
om 2@ j 
+ (8 cos? Он — 1) c ae > (309 qf? — »»| 8i... (4.9) 


where 000 i is the cosine of the angle formed by axes X® (of the system attached to 
the path И. and axes xe i that attached to the q-th molecule). 
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On effecting the transition from (4.4) and (4.5) to (4.8) and (4.9) we had assumed 
a Lorentz field as appropriate for the light wave E? within the sphere S, i. e. 
24.2 
ee E Е. 


The Kerr constant is defined as 


fremit СОУ a ee mn] 
D n E? бп? п +2 № +2] Е?’ 


where n; and n, are the optical refractive indices for a light wave field Е“ parallel 


(4.10) 


and perpendicular to the electric vector E. 
From (4.8) and (4.10) we obtain the following expression for the Kerr constant 
of a dense medium: 


га NR NS Ns п? 4-2 : Es 2 oe oe 
= ТЕ Vs (^ 37 | Е D Зеу — Сил 


і ij 


+ 2H (зыў) Qi og + < Qi (Goff of? — 1) > + 


Ns 
n: Y qr 
Tg pi 3 Y (Bae? ee o% a (4.11) 


q=1 r= 


Analogously to the Kerr constant, we shall define a Cotton-Mouton constant by 


т=п 1 S cese Cd 1 


С" = 


n H? 6n? nj-2 т +2 Н? 


and eq. (4.9) yields the Cotton-Mouton constant of a dense medium: 
4 л № {п? +2 A ede о о 
C = is Y. 2) Ч 2. Seg — Gig + 
ed айай асу (3 O рУ) [ M (4.12) 


Neglecting T oc. we obtain from (4.11) and (4.12) the following relations 


holding for molecules presenting axial symmetry: 


n +2 oe oe 
K* әт TRE 3n ) (Es 2 e Rom + OF Re); 


and 


А п + 2 от 4.13 
= و‎ | 3л } (48) ө Rem; (4.13) 
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where 


oe 2 е е о o 2 5 С 2 
Giga 45 kT (азз — а11) (азз — ап); 92 = 15 PT? (азз — аз) иѕ, 


2 т т о 0 
@ = 15 LT (a33 — алл) (азз — a11), | (4.14) 
and the correlation factors: 
Ns T 
1 
Rom = 3 d ; (3 cos? Ө, == 
-— 
Ns №5 
1 ps 
Ек = o: < 124 (3 cos O,, cos Op, — cos Өл.) 7 (4.15) 
q= r= 


with @,, representing the angle formed by the axis of symmetry of molecule p and 
that of molecule q. 
On assuming A. Piekara's model discussed in 83, eqs. (4.15) yield 


Rem =1+5 COR = EE (1—22), 


Кк =1 +3 (coos 95 —1) +2 (cos 04,» = 1 + ( +32 +21 (4.16) 


which are A. Piekara's (1939, 1950) correlation factors appearing іп the Kerr and 
Cotton-Mouton effect theory. 


5. On the dielectric constant of a diamagnetic dense medium deformable in electric 
and magnetic fields 


The total potential energy of a system subject to the effect of a strong electric E 
and magnetic field H forming the angle Q is (apart from mechanical deformations): 


U(r, Es, Hs) = U (т, O) — M, Е, - AEE TAHT 
1 ee 1 е 
єз! (Bi E, Е, Е, +3 Ве", Es HH, + ...) — (5.1) 


X 1 (Core Е, Е, E, E, is бот Ec Er H, Но zia E == ... 3 


where U (т, 0) represents the internal potential energy of the system, М, — the electric 
. moment of the macroscopic sphere S, 4%, and A™ — the tensors of the electric and 


diamagnetic polarizability of the macroscopic sphere 5 respectively, В“, С“, — 
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the electric non linear deformation tensors of the macroscopic sphere S and, finally, 
Be» and C77 „— the electro-magnetic non linear deformation tensor. MPA ASB, 
and Covo are tensors symmetrical in all indices, Bey is symmetrical in т and >, 
whilst Cv, is symmetrical in the separate pairs o, т and v, 0. As above, (5.1) 


and (2.3) yield 


бот е М.М. т 
2 (с, Es, Hg = Са T | (ж. i мм. ES rin 


- L 


uU 3۰5! ET (до: д, ^r до» Ore ME доо бу») lee zi kT (342, А» T 4M, BA 2E 


М.М, е M, M, 1 
uu (в, X TE J] ES + ТЕТ @— cos Q) er Öne + 


1 e m т 
+ (3 cos? 2 — 1) (05, б + до 0,,)] х [= лг LT (Aor А» + 2M, Brive) + 


М, М, А" 
ae d ES H$ + © : (5.2) 


By (5.2), eqs. (2.1) and (2.2) yield the dielectric permittivity of a dense medium in ge- 
neral form: 


JE S Cans 
MES (s nate (s Mt) ву) + 
p as cms + 2С, + ET [47M 7 (Bo SP RUE АА Е 
zb um | 2 Ао М, М, + 4 45, М, M, + уы M,M,MzM,) > ae 
+28 Co cost 0) [ee ema pest (4: + Sew + 


em е с: : 
a (3 cos? Q — 1) fee "T Tm M, Bio: = ET (4 "js аң > xk a (5.3) 


Neglecting the deformation terms with polar molecules presenting axial symmetry 
(5.3) yields: 


RE Ns 9 Es us 3 cos? Q —1 O” R O” R ) H$ + 
ع‎ l= 2 Fs LEL ( 1 см + 03 LK, 


+ (Of Rem + 207 Кк — OF Rs) #8 , E (5.4) 
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where 
2 e е т т ет 2 т 2 
9i = ET (азз — ал) (a3 — аз); 92 = 45 PT? (a33 — a3) из, 
Ө x. 2 (аз m ot )?; Ө“ m 2 (a3 ауу) us z OF == ЗЕ Е (5.4а) 
о З 45 k3 T3 


and a, — З (2 af, + ag) represents the mean electric polarizability of а molecule 
within the sphere S, the correlation factors Roy and Ry being given by (4.15) and 
Rp and Rg by (3.18). 

For the variation of the dielectric constant in a strong electric field, eq. (5.4) 


yields 
Aci, = 120 75 (Of Rom + 20$ Вк — OF Rs) dle) ЕЗ, (5.5) 
: | 


where d(e) is defined by eq. (3.92). 
Similarly, for the variation of the dielectric constant in a strong magnetic field, eq. 


(5.4) yields: 


Meee, = (3 cos? О — 1) (6*" Rom + Ө Rx) d(e) HÀ. (5.6) 
S 


The microscopic method yields identical results for the correlation factors Rp, 


Roy, Ёк and К; (A. Piekara and S. Kielich, 1957, 1958). 


The author wishes to express his most sincere gratitude to Professor A. Piekara, 
Head of the Laboratory of Dielectrics of the Institute of Physics of the Polish Academy 
of Sciences in Poznan, for suggesting the present subject and for his valuable advice. 


КРАТКОЕ СОДЕРЖАНИЕ 


. С. Келих. Семимакроскоповое представление теории нелинейных явлений в диэлектри- 
ческих жидкостях подвергнутых сильным электрическим и магнитным полям. 


Автором построена общая полу-макроскопическая теория нелинейных эффек- 
тов молекулярной ориентации, вызываемых в жидкостях сильными электриче- 
скими и магнитными полями. Автор дает вычисление отклонения от квадрати- 
ческого эффекта диэлектрического насыщения в сильном электрическом поле 
для полярных жидкостей с изотропно поляризуемыми молекулами ($ 3). Затем 
дается вывод общих выражений определяющих влияние сильных электриче- 
ских и магнитных полей на показатель преломления ($ 4) и диэлектрическую 
проницаемость ($ 5); основываясь на этих выражениях, автор вычисляет по- 
стоянные Керра и Коттона—Мутона, а также изменение диэлектрической прони- 
цаемости в электрическом и магнитном поле. Молекулярная интерпретация этих 
результатов приводит к формулам полученным раньше молекулярным методом 
(см. А. Пекара и С. Келих). 
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THERMOEMISSION OF ALLOYS* 
By КО, Е. WOJCIECHOWSKI 
Department of Experimental Physics, Wroclaw University 
( Received November 19, 1957) 


A semi-phenomenological theory is presented for thermoemission from binary 
and ternary alloys which form a superlattices at low temperatures. It has been found that 
the thermoemission current has a local anomaly in the neighbourhood of the order-disorder 
transformation temperature. The character of this anomaly (maximum or minimum) 
depends on whether the order-disorder transformation is a first-order or second-order 
phase transition. Richardson's formula has been derived for the alloys undergoing an 
ordering process; Richardson's constant 4 for the case of these alloys is a function of the 
parameters of long-range order. It is noted that the thermoemission current for alloys 
with a non-cubic crystal lattice is anisotropic. 


1. Introduction 


It is known that many physical quantities such as electrical resistivity, specific 
heat, diffusion coefficient, and others have, in the case of binary alloys with a super- 
lattice at low temperatures, a local anomoly in the neighbourhood of the order-disorder 
transformation temperature Т. This anomaly is caused by the effects of the ordering 
of the atoms in the crystal lattice. In the case of the specific heat (e.g. Sykes and Jones 
1938), of the electrical resistivity (e.g. Sykes and Evans 1936), or the diffusion coefficient 
(Krivoglaz and Smirnov 1955), or many other physical quantities these anomolies were 

adequately investigated from the experimental and theoretical sides (e.g. Fowler and 
Guggenheim 1939). 

It may be expected that the thermoemission current too will show, at the point 
order-disorder transformation, an anomoly caused by disappearance of long-range 
order in the alloy. 

In thc present paper, a semi-phenomenological theory of ATUM AES from 
binary and ternary alloys is given; the variation of the thermoemission current as 
a function of the degree of long-range order is also discussed. 


* This paper was read at the 15th Congress of Polish: Physicists in Wrocław on November 
11221057 © . 9 А 
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2. Binary alloys 


In order to calculate the thermoemission currents we make use of the well-known 
expression for the density of the emission current from a unit area of the metal surface: 


m Ux dv, 
es: Bl fo је ka [E — КУШ! ©) 


where о, Vy, V, are the vector components of the electron velocity, Ep the Fermi 


energy, E the energy of an electron moving through the crystal lattice of the alloy. 

Smirnov (1947) found the relation between E and the parameter of long-range 
order. For the sake of simplicity, we shall examine alloys having a body-centred cubic 
lattice. For such a lattice 


7 
Е= В + [a — с)? ys? + 640 cos? um cos? p: cos? Зр (2) 


where B, y, q, £g are constants, c the concentration of one of the components of the 
alloy, a the crystal lattice constant, s the parameter of long-range order. 
Eq. (1) can be rewritten in the form 


e 2eh ky dk, 


dk, 


+оо 
dk, 


h 
wheie m* = Jap is the effective electronic mass, which may be calculated by means of 


E 


(2); k,, ky, k, are the components of the wave vector К; and k? = ls uot ig ig. la). 


Performing the integrátion in Eq. (3) in the same way as for the derivation of 
Richardson's classical formula (see e.g. (Kaptsov 1955)), we finally obtain 


Lo 4mekim*:,,, _ Ра ВЕ 
Ps a De us TANE (4) 
The above formula for the density of the thermoemission current from a binary 
alloy agrees in form with Richardson's formula, while it differs from the classical 
Richardson formula only in that 4 and W, depend on the long-range parameter. 
Since 
| Agtek®m* 


A t в 


Е 9E (3) 
ep | ak К — : E) | «| +1 : 


"y #2 
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then, calculating from Eq. (2) the effective mass for the upper, nearly filled band, 
we obtain 


a l(q — ©)242s? + 6420] ® (5) 
50а“ 


W , may be obtained as a function of s in straightforward way by expanding W into 
a series in 5 in the neighbourhood of the point of the order-disorder transformation, 
that is, in the neighbourhood of the point s — 0: 


; IW 1 (921. 
Р CETERI | e а 2 
a(S) = W,(s)o 4 | F К 5 | I E 2 зо (6) 
Since the linear term vanishes because of the symmetry, we obtain, by considering 


only the term containing s?, 


W, (5) = W, (o) + as? (7) 


a 


where x = ые 
: ET D 952 m 


Denoting the work function by ф we obtain on the basis of the above considerations 
p(s) = W,— E, = W,(o) — Ep + «s? (8) 


A similar relation was introduced in the paper of Sokolov (1951), but the constant 
appearing there differs from the above ones, since Sokolov assumed that Ep too is 


a function of s and he expanded the difference W (s) — Eg (s) into a series in s making 
use of the well-known formula! E, = 


Since this assumption, however, is incorrect (5 = 1 for T— 0), Sokolov's formula 
only accidentally has a form similar (as regards the constants) to the correct formula (8). 
It should also be noted that Sokolov assumed, on the basis of the general theory of 


ҮШ IW ? i 
second-order phase transitions, that 5 i) <0 below the Curie point. It is known, 
5 0 


however (Lifshits 1941), that the order-disorder transformations are not always second- 
` -order phase transitions, sometimes they are first-order phase transitions, and then 
the temperature of the order-disorder transformation does not coincide with the 


927. А 
Curie temperature; hence the case 5 ;]29 should also be examined. On the 
s" Jo 


basis of the above we may finally write Eq. (8) as follows: 
g(s) = (o) ± a -— (9) 


where 


p(o) = W, (0) — Ep _ Q0 


1 Private communication from A. V. Sokolov. 
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denotes the work function at the point of order-disorder transformation (s — 0) 
(Fig. 1а). 
On the basis of (5) and (9) we can rewrite (4) in following way 


(s) 


9 
I(s) = q — cs + 64eg] T?e вт (11) 


Aha} al 

It may be seen from (11) that if the order- disorder transformation is a first-order 
phase transition that is, if T, 7^ T, (T, is the Curie temperature) then the thermoemis- 
sion current has at T, a fet maximum (Fig. 1c) and if the order-disorder transforma- 
tion is a second- iue phase transition that is if T, = Т, then the thermoemission 
current has a local minimum at T, (Fig. 1b). We then see ‘that the effects of the long- 


J(T) 


Was) 


ЕНСЕ 


= 
27 fo) Fle f€) 
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Fig. 1. 


-range order parameter s in the factor before the exponential expression and in exponent 
of Eq. (11) are superimposed on one another and produce an anomaly stemming from 
the ordering.” : 

It is true that Eq. (11) was derived only. for alloys having a body-centred zi 
lattice, but we shall show in Section 4 that it can also describe, to a sufficient approxima- 
tion; the influence of quss on the thermoemission from alloys with any type of crystal 
lattice. 

In order to investigate the effect of order on the phenomenon of thermoemission 
the temperature in the neighbourhood of T, should be examined. A first suggestion 


` 2 Local discontinuities of the thermoemission current curves at the melting points of metals are 
known (e.g. Goetz 1927). 

` 3 Parameter s is a function of the temperature (Guggenheim and Fowler 1939), and therefore in - 
order to obtain the above results it may be expressed by the temperature in the рынке of T, and 
inserted into Eq. (11). 
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which comes to the fore is that any eventual influence of this effect on thermoemission 
can be taken into account only for those alloys for which T, is sufficiently high; for 
alloys characterized by a low transformation temperature; there is no sense of investiga - 
ting the way in which the order produced by an increase in temperature influences the 
thermoemission, because the latter is negligible at low temperatures. 

There exists, however, a large number of binary alloys having a rather high trans- 
formation temperature; we mention here some of them: CoFe — Г = 800°C; CoPt — 
T, = 825°C; CrFe — Т, = 815°С; CuPt — T, = 812°C (Oriani 1954); Fe;Pd, — 
— T, = 780°C; FePd, — T, = 780°C (Hultgren and Zapffe 1939). Investigation of 
the effect of order on thermoemission in the case of the above-mentioned alloys would 
be fully justified. 


3. Ternary alloys 


We shall now consider thermoemission from ternary alloys. From the mathemati- 
cal point of view this case is more complicated, but from the physical point of view, 
as we shall show, it is analogous to the previous case. : 

In order to describe long-range order in a ternary alloy four parameters of long- 
-range order: sj, Sg, 53, s, should be introduced (Wojciechowski 1956); thus the energy 
of an electron moving through the crystal lattice of such an alloy is a function of these 
four parameters: E = E (51,55,53, 54, К). Making calculations similar to those for 
binary alloys, we find that in this case Richardson's constant A depends on four 
parameters of long- range order A = A (51, 5, 5з, 51): Since, however, chiefly the expo- | 
nential factor influences the thermoemission current anomoly, if we proceed as in the 
case of binary alloys, we may derive a formula giving the work function in terms of 
the parameters of long-range order: 


4 


1 [89*W, 1 227, 
Ф(51, Sos $3, $4) = Ф(0) + 2 э | ER } = 5 (5528), 5 5 -h (12) 


i=1 


Since only the neighbourhood of T, interests us, we can assume that 51 ~ sg ~ 


~ 53 ~ 54 = 5, from which, according to whether T, + T, or T, = T, we have 


Р) = ро + а, | | (18) 
where 
А | 
T 92g =- 927, 2 
“=> | 24 \ 4 Ds ma |. 
{= Ч i*#j 


We have thus obtained the same result as in the case of binary alloys. It may 
therefore be assumed that in the case of ternary alloys the thermoemission curves іп 
the neighbourhood of T, has the same shape as in the case of binary alloys. 

D 
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4. Significance of the one-electron gas model and the method of effective mass for 
investigating the influence of order on thermoemission 


The classical Richardson formula for the thermoemission current density for 
pure metals was derived on the basis of the Sommerfeld model of a free electron gas. 
In this formula the constant 4 is a universal constant. All the above considerations 
concerning thermoemission of alloys are based on the one-electron model and make 
use of the effective mass method. 

As we have seen, it follows from these considerations that A is not a universal 
constant, but a quantity dependent on the long-range order parameter s and on the 
concentration c. This difference is based on the fact that in the case of pure metals it 
may be assumed as a first approximation that the electrons move in a constant potential 
independently of the temperature, while in the case of alloys such an assumption is 
completely without sense. Electrons collectively moving through the crystal lattice 
of the alloy are continually in a field of variable potential, which is particularly impor- 
tant when the atoms of the alloy spontaneously change their position and form a disorde- 
red structure near the order-disorder transition temperature. In order to take into 
consideration the interaction of the electrons and the crystal lattice of the alloy, we 
have introduced an effective mass, which in the first approximation allows the interac- 
tion to be taken into account. 

As it is know, the effective mass is, in general, a tensor with components. 


O2F\-1 / 02E \-1/ O2E үа 
Ej Е в 
2p \-1/2р\-1 ER ЛЕТ 
TERM TEC Eom ой (14) 
9k, 3 ky aky Ok9k,] -— 
ОЕ \-1/ д2р \-1 [д2 \=1 
ea ne (S 
and therefore depends on the crystallographic direction for all structures except cubic 
ones. In the case of a cubic structure all the components (14) are equal. For other 
structures, however, the entire matrix (14) should be considered, which, as a result, 
would give three different effective masses along the x, y, z axes; chosen in the direc- 
tions of the principal axes of the effective mass tensor: m;, т», m;. In order to са]- 
culate them for alloys with non-cubic lattices, use should be made of the general 
expression (15) from Smirnow’s paper, from which it follows that the effective mass m7 
will depend on the paramater s in a different way. This only reflects the dependence 
of Richardson's constant on s, while it does not at all influence the exponential factor. 
Since it is primarily this factor, however, that causes the jump in the thermoemission 


current in the neighbourhood T, the character of the curves in Figs. lc and la for 
any alloy remains the same. j 
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Also noteworthy is the fact that on the basis of the above considerations concerning 
the effective mass, it should be expected that for both binary and ternary alloys with 
non-cubic lattices, we will be dealing with anisotropy of the thermoemission current. 
This anisotropy reflects only the size of the jump in the thermoemission current at 
the transformation point, but does not, however, change the character of the current 
curves. Thus, although we made detailed calculations only for alloys with cubic lattices, 
the results obtained, while applying only to the shape of the curves in Figs. 1c and 1b, 
have a general character and apply to alloys of any composition and with any crystal 
lattice. 


5. Conclusion 


The above results make no claim to be a qualitative solution to the problem of 
influence of the ordering phenomenon on thermoemission of metallic alloys and were 
derived. on the basis of an approximate physical model without confrontation with 
experiment. It therefore seems worthwhile to make an experimental quantitative check 
of the results of this paper. 

The thermoemission current anomoly at the order-disorder transformation point 
for binary alloys was discovered independently of the present author by A. V. Sokolov; 
the results of his work, however, have not been published (private communication). 
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КРАТКОЕ СОДЕРЖАНИЕ 


К. Войтеховский, Термоэмиссия сплавов. 


В статии дана полуфеноменологическая теория термоэлектронной эмисии из 
двойных и тройных упорядочивающихся сплавов. Найдено, что ток термоэлек- 
тронной эмисии имеет локальную аномалью в окресности точки перехода порядок- 
-беспорядок. Характер этой аномалии (максимум или минимум) зависит от того 
является ли переход порядок-беспорядок фазовым переходам: І или П рода. 
Выведена формула Ричардсона для упорядочиващюихся сплавов. Найдено, что 
постоянная Ричардсона А является для этих сплавов функцией параметров 
дальнего порядка. Обращается также внимание на анизотропность тока термо- 
электронной эмисии в случае сплавов имеющих некубические кристаллические 
решетки. : 
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INVESTIGATION OF THE EFFECTIVE MASS OF CURRENT 
CARRIERS IN GaSb 


Ву 7. КОРЕС 
Institute of Physics, Polish Academy-of Sciences, Warsaw. 
( Received December 30, 1958) 


The thermoelectric force, Hall effect, and electric conductivity for-3 samples of 
GaSb-p and 1 sample of GaSb-n were measured in the temperature region 200?—400*K. 
The scattering of carriers on phonons and ions were taken into account in calculating the 
values of the effective mass, the latter proving to be a function of both the concentration 
and the temperature. 

Additional experiments on the measurement of the thermoelectric force and heat 
conductivity for one sample of the p-type and one of the n-type were made at lower tempe- 
ratures (attaining the temperature of liquid helium) indicate that for large current carrier 
concentrations the values obtained for the effective mass at the higher temperatures are 
not affected by the phonon drag effect. 


By taking into account the mechanism for the scattering of free current carriers, 
the effective mass of electrons or holes can be determined from measurements of the 
conductivity and thermoelectric force in semi-conductors. 

* The mechanism of scattering in semi-conductors with a structure and bonds of 
the germanium type is relatively the best known. For valence semi-conductors there 
exists, in general, a temperature region in which the dominating role is played by scatte- 
ring on phonons and ions. The values of the effective mass determined in this region 
proved to be in quantitative agreement with the values obtained by other methods, 
in particular, on the basis of the cyclotron effect (Kipp and Kittel 1955). Analysis of 
these values for Ge, and more distinctly for InSb-n suggests that the effective mass 
of the carriers changes in a regular way with changes in their concentration and tempe- 
rature, the effective mass increasing with an increase in temperature and increasing 
even more so with an increase in the carrier concentration. 

In order to compare the above-mentioned materials with other, similar ones, 
a series of measurements were made of conductivity, the Hall effect, and-thermoelectric 
force on gallium antimonide GaSb. With the use of the method described in a previous 
paper by the author (Kopeć 1956), a group of these values permits the determination 
of the effective mass of the carriers. The energetic gap between the fundamental band 
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and the conductivity band is relatively large in gallium antimonide and equals AE 
— 0.78 eV (T = 300°K); we may therefore use here the simple expression for thermo- 
electric force a for the non-degenerate case 


d Hes [Se wi | , 2@л&Тт*)у' | Е 
vi em lá E du nh? | (1) 


is the Boltzman constant, 


is the electronic charge, 
is a constant dependent on the scattering mechanism, 


Soho m 
Ф 


is the carrier concentration, 
olo; is the ratio of the real resistance to the resistance with phonon scattering. 


If we take into account the fact that at the temperatures Т = 200°К and Т = 
= 400°K the carriers are scattered chiefly on phonons and on ions, and if we employ 
the values A (0/01) and r(o/o;) calculated by Anselm and Klatshkin (1952), and formula: 
R = r(o[oz) Те, we find the following values of m* for different samples of GaSb. 


Table 1 


Values of effective mass in electronic mass units 


Samples of p-type 


0.39. 0.19 0.15 


Sample of n-type 
T(°K) 
п 12 x 10:3 сш 


200 


Fig. 1 shows how the thermoelectric force «, Hall constant R, and the electric 
conductivity с depend on the temperature for one of the GaSb-p samples at tempera- 
tures of from 200°K to 400°K. — i | 

Table 1 indicates that in contrast to germanium and indium antimonide the 
effective mass of the holes in GaSb is greater in pure samples than in the more impure 
ones, and therefore has a greater concentration of carriers. 

This fact — unexpected on the basis of the purely band model — suggests that 
the scattering mechanism has not been properly taken into account. 

As may be seen from Table 1, a particularly large jump in the value of the effective 
mass is observed in passing from a sample of average purity to a pure sample. This 
fact suggests that a role may be played here by the “phonon drag” effect which was 
not taken into account in formula (1) and connected with the interaction of phonons | 
not being in equilibrum (i. e. of mean velocity in the direction of the temperature 
gradient) with the electrons or holes of the semi-conductor. (Gurevich 1945, Frede- 
rikse 1953). It is now already known from numerous semi-quantitative theories of 
the phonon drag effect that the effect causes an additional increase in the thermo- 
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electric force a. This effect is great for very pure samples, and disappears for samples 
of considerable impurity. In pure samples it causes the appearance of a resonance 
maximum at low temperatures. As the carrier concentration increases the maximum 
spreads and shifts to the higher temperatures. It is therefore not excluded that the 
anomoly in the large value of m* for purer sample is connected with the fact that the 
phonon drag which already begins at these concentrations was not taken into account. 

Further investigation of the scattering mechanism required that the experiments 
be carried out at low temperatures. For this purpose a special experimental set-up 


cm3.coul 


100 


ج 1000 
2 4 3 2 1 


Fig. 1. Thermoelectric force a, Hall effect R and electric conductivity o for a p-type samples 


was constructed for measuring the thermoelectric force and' for the thermal 
conductivity. 

In this set-up measurements were made on the GaSb samples at our disposal, 
one of the n-type and another of the p-type. | | 

Fig. 2 shows the general scheme of the low-temperature experiments. The sample 
was placed in a copper thermostat from which the gas was pumped out to a pressure 
below 10-4 mm Hg. The thermostat was immersed in a vessel containing activated 
carbon with a well-extended surface. 

By pumping out the helium adsorbed on the carbon it was possible to obtain 
any constant temperature in the range of from 5°K to 70°K in the course of an hour. 
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This simple arrangement permitted the taking of measurements under stationary 
conditions over a wide range of low temperatures with the use of only liquid helium 
and liquid nitrogen. The temperature of the ends of the sample mounted as in Fig. 2 
was measured for temperatures below 25°K by means of carbon resistance thermome- 
ters whose sensitive elements were plates about 0.3 mm thick cut from carbon mass 


Distribution 
table 


Liquid. helium 


Jessel with 
activated carbon 


Fig. 2. Set-upf or measurement of the thermoelectric force and thermal conductivity at low temperatures 
resistors of the radio.type. The ends of the sample were set in a copper block inside 
which were the thermometers and their leads. 

This method of mounting the thermometers guaranteed the accuracy of the tempe- 
rature measurements. The thermometers used were unusually sensitive, the resistance 
of such a thermometer at a temperature of 300°K being about 100 2 at 77°К — 
140 Q, at 4°K — 8000 О, at 1.5°K — 80000 О. Above 25°K it proved more convenient 
to use copper-constantan thermocouples. 

Fig. 4 represents the results of measurements of the thermoelectric force for 
samples of the p-type and n-type and containing large amounts of impurities. The 
thermoelectric force increases monotonically with an increase in temperature for both 
these samples. For samples of the n-type the curve of the variation of the thermoelectric 
force with temperature is similar to that in the case of metals; we are thus dealing 
here with a case of carrier degeneration. On the curve representing variation of the 
thermoelectric force with temperature we observe a discontinuity in the vicinity of 
100°K. Without additional measurements for samples of various concentrations of 
free carriers it cannot be determined whether this discontinuity is evidence only of 
a transition to a region of degeneration or also indicates a.drag effect. It can only be 
stated that even if this effect occurs here, then it is slight and does not disturb the 
monotonic increase of the thermoelectric force with temperature. 
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r Suspension 


Termocouples 
Thermometer leads 


Copper housing of the 
thermometer 


Sample mounting base 


400 


200 


200 


i Fig. 4. Thermoeleçtric force a at low temperatures for p- and n-type samples 
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Fig. 5 shows the results of the measurements of thermal conductivity, which for 
the GaSb-n and GaSb-p samples agree within the limits of experimental error (10%). 
This fact indicates a uniform and probably not great contribution of electrons to the 
heat conductivity and confirms the conclusion concerning the insignificant role of 
the drag effect which should be reflected in the curves of the thermoelectric force 
and thermal conductivity (and therefore in other quantities connected with transport). 


Мест ies sper 
1.5 
1 
0.5 
- 100 200 300 400 T*K 


Fig. 5. Thermal conductivity x at low temperatures for p-type samples 


The thermalconductivity curve at temperatures above the maximum indicates that 
phonons are scattered on phonons; at the lowest temperatures measured, the phonons 
are scattered on the boundaries of the samples and on the impurities. 

The presented data indicates that for GaSb samples with large impurities the pho- 
non drag effect is slight and therefore does not distort the values obtained for. We 
have the hope that further measurements made on pure samples will allow the investi- 
gation of this effect on the measured quantities to a greater extent. 

An interesting result of the investigations is the establishment of a considerable 
increase in the effective mass with an increase in temperature. The increase in the 
energy of the carriers, which with an increase of temperature pass into the interior 
of the band could only cause an effect constituting 25% of the observed effect (as 
follows from an estimate assuming even a considerable change in the energetic structure 
of the band with the filling of states deeper and deeper into the Band such as in InSb). 
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It should therefore be supposed that the increase in temperature causes a considerable 
modification of the same band. 

I would like to thank Professor A. F. Joffe for his interest, valuable discussions, 
and aid during this work. 

Part of the work connected with the measurements at temperatures below 77°K 
was made in the liquid helium laboratorium of the Physical-Technical Institute in 
Leningrad. I would like to thank the administration of the institute and the laboratory 
chief N. M. Rejnov for making my work possible and for his effective aid and I. Tim- 
chenko for his help in the measurements. 


КРАТКОЕ СОДЕРЖАНИЕ 


3. Копеть, Об исследовании эффективной массы носителей тока в GaSb. 


В области температур 2009 К — 400° К измерено термоэлектрическую силу, 
эффект Холля и электрическую проводимость для трех проб GaSb-p и пробы 
GaSb-n. Учитывая рассеивание носителей Ha фононах и ионах подсчитано зна- 
чение эффективной массы являющейся функцией концентрации и температуры. 

Добавочные опыты над одной пробой типа р и одной типа п касающейся 
измерения термоэлектрической силы и проводимости теплоты сделанные в более 
низких температурах (достигающих гелевым) указывают, что дяя больших кон- 
центрации носителей тока полученные в температурах высших значения, эффек- 
тивных Macc не нарушаются эффектом фононового увлечения. 
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A SIMPLIFIED METHOD OF SOLUTION FOR A RANDOM 
WALK PROBLEM OF A FEW UNEQUAL STEPS 


By Denis O’ CONNOR 
Institute of Nuclear Research, Warsaw ` 
( Received January 27, 1958) 


A simplified method is presented of calculating the probability that the amplitude 
of the resultant vector formed by the vector addition of a finite number of vectors of known 
amplitude and random orientation, will lie in a given range. This method is simpler to 
use for small numbers of vectors unequal in length, than methods known from the 
literature. Recursion formulae are given for calculating the probability functions for 
assemblages of vectors of equal length. 


The problem we are concerned with here can be stated as follows: N ‘vectors of 
known length and of isotropic distribution add to give a resultant R. We wish to find 
the probability P(R)dR that R lies between R and В + dR. This is equivalent to 
finding the probability that a "walk" in space of N steps of known length and random 
mutual orientation, will end at a point between R and R + dR distant from the origin. 
By making appropriate assumptions about the amplitude distribution of the vectors 
this problem can be applied to, for example, Brownian motion, calculation of the 
resultant force on stellar bodies etc. (Chandrasekhar 1943) More recently it has been 
applied to the calculation of the probability that a nucleus, excited in a (my) reaction, 
will escape from the molecule in which it is bound as a result of the momentüm impar- 
ted to it by the photons emitted (Cobble, Boyd 1952, McCallum, Maddock 1953). 

Rayleigh (1919) has given a solution to this problem which has been elaborated 
by Chandrasekhar, and which can be arrived at by ee use of the unique pro- 
perties of the Dirichlet integral. The general solution is: 


рк) === i f. sin (|pl: вр {fF del 2 on (8. -D paie | 


]== 1 


Where J, is the length of the j-th, vector and R is the resultant of the addition 


of the N vectors. Р 
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It will be appreciated that the evaluation of this integral is very tedious for even 
a small number of vectors, if all the vectors J, are of different amplitude. This will 
be the case in the (ny) reaction problem mentioned above. In the following a method 
of solution is given which, while leading to the same results uses simpler mathematics 
and is more easily applied to the case of a few vectors of different amplitude. 

We first consider the problem of two vectors A, and kọ which cómbine to give 


a resultant r: 


r= hy rh 
r? = k? + k? + 2k cos Ө 
whence rdr = — kk, sin OdO (2) 


Now $sin O dO isthe probability that k, will make an angle between Ө and 9 + dO 


with hy. It follows from (2) therefore that the probability P(r)dr that r lies in the interval 
r, r 4- dr is given by: 


үз 
= — —d 3 
P(r)dr ge (3) 

The basis of the present method is the repeated use of equation (3) by breaking 
an JV-vector problem down into a series of (N—2) 2-vector problems. 

Suppose P,(r,)dr, represents the probability that the amplitude of the vector 
resulting from the vector addition of vectors k T a lies between r, and г, + dry. 
We shall show in the following that there exist a number of different expressions for 
P,(r,) which are applicable in different ranges of r,. Suppose that in the range of г, 
from r,, to г the appropriate expression is P7(r,). Nav suppose we add another vector 
k, +1 to 7, to obtain г, , ү. If r, were constant we should have according to equation (3): 

T 
Taas aN Ta p 


n-4- 1n 


To take account of the probability distribution of r, we must modify this expres- 
sion to: : 


rb 
Т. P; (Tn) 
j2. cec agi. п\п) 
+1 (+1) Tha ff mae ; (4) 
"а 
where the integration is taken over all values of г, which when combined with A, |, 


will give the required values of r, |. Now the maximum and minimum values of r, for 
which the expression P, (r,) is valid are г» and rj respectively. We also ed 
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It follows that 


From the inequalities (5) it follows that the values of r, and г, to be used in (4) 
depend on the range of values of r, ,, being considered. For example if r, |, 4,4, is 
greater than г,» the maximum value of г, will be T, = Г. Reasoning in this fashion 
we obtain: 


Tul aa» Tact бы УЕ р 
new too tni cp AINT 
Ла fip Et xy ees 
EE | Jm] 2 kb, a |, | E TT Eia | = Tal 


which can be rewritten in the following form: 


Tati Ао ва [AES Ty = Ги 
Tata < Га — Ёр и 
И 
ELEM (6) 
Falc ka Regn 
ntl 2 Tag i LET a res LA 
Tapi E LET EU Tak LET ES 


In addition we note that the upper limit on г, ,, is г, + А, у and the lowest value 


assumed by г, ұу will be r,, — k,,, for Г > А, and zero for г < Ё, 


п+1 
Thus for гд > Ё, 
| La uel a Tat Kets ба) 
and for гд < Ё, | 
га + nyt > Гл > 0 7 (7b) 


The value of P, , ,(r, ,,) in the range of values of г, ,, given by either (7a) or (7b) 

will be found by substituting in equation (4) for r, and r, from the relations (6). From (6) 
and (7) it can be shown that there will be three or four different sets of limits r,, ть 
depending on the relative values of Ги» Тв and Ё, ул. It follows that from the single 
expression P’(r,) there result three or four different expressions for P, +1 (71) which 
are applicable in the different ranges of values of г, , resulting from the consideration 

. of the inequalities (6). The assumption made above with respect їо P,(r,) is thus justi- 
fied. It will be noticed that the ranges of г, үү are formed by adding k,,,, to, or sub- 
tracting it from the limits of the ranges of r, : 7,1, Т. By induction it follows that т 
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where @ is a positive whole number even or odd depending on j and m. 

To evaluate P_(r,) using this method it is obviously necessary to have evaluated 
P, ar, a). P, sr, 9) ... P2(T2) or in other words to start from P (72) which is known 
(equation 2) and work upwards evaluating all the P(r) for all their subranges of r. This 
is in effect what is done in integrating expression (1) by parts, but the present method 
is simpler in that the expressions to be integrated are simple algebraic expressions 
in powers of r, the crux of the method being the determination of the limits of integra- 
tion by the use of the inequalities of the for of (5). It is not necessary to find explicitly 
the а in equation (8) since the values of г, automatically result from consideration 
of the inequalities used in determining the limits. As an ilustration of the form taken 
by the P(r) we give below the solutions for a three vector problem, where we assume 
leg М. 


0 — ry < № + А — А, Р-з(@з) = Е 
r«(r. ky th + Е 
hy + h-hench hk псы E 


kı + kg — № < в < В + kg — №» P3(73) = Е 
Tg ( + ka + ks — гз) 


ky + № — Ез < r3 < kı + ka +, Pars) = Ak kk, 


We have here also made the assumption that 0 < ka + ka — А < kı +h, — 
— ka < k, + ka — Ез, which will not always be the case. In fact appropriate solutions 
can be given only when the relative values of kı, А, and kg are specified. 

For the case where all vectors are of equal length, k, some simpler expressions 
can be given. First we notice that the limits to the ranges of r, r,, are all multiples 
of k. Equation (8) becomes: 


Now if we have an even number of vectors 2n, the smallest non-zero positive 
value for г, will be given when in the summation 2(—1)* there are n + 1 even values 
of @ and n — 1 odd values, in which case Г = 2k. The next value of r,, Will be for 
п + 2 even and n — 2 odd values of a i.e. г, = 4k. Thus for 2n vectors it is seen that 
the r, are given by 


2k, 4k, 6k etc. 
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and for an odd number of vectors it can be shown that the r,, are given by: 
ks, OR SE sete. 


That the fundamental interval between the Ty, must be 2k obviously results from 
the fact that a single vector Ё makes a contribution to the resultant over a range of 2k. 


Let „Ро„(т„) be the expression for the m-th range of ra, (counting from zero) 
for which we have: 


2k(m — 1) « ry, < 2km 


Now rs, can be derived from 7,,,_ by adding a vector k. It follows that Ps, (rs, 4) 
makes a contribution to , P5, (r,) for a range of rj, , given by: 


3 1 
2k (m — z) I 24 (m -- ы 


In this range of 7, , there exist two expressions for P$, j(r,, ,) i.e. 
3 1 
me en—\(Toen24). 2b | m — э] < Tan- < 2k | т — P 


1 1 
m+1Pon—y (Tan-1), 2k (m Ж 3 < Tan-1 < 2k (m ie 3 


We can therefore write, using equation (4): 


Tb, Tb, 
_ Ten т+1 Pon-1 Ж) Ton mPon—1 (ran—1) d 10 
mP an (Ton) = әр Saat peepee, а-у + әр аге жуа Сов (10) 
та Таз 


where we must now determine the limits. From the conditions on г, given above 
we have: 


1 
Tan +k 2 ты < 2 (m +2) 


1 
(ran — k| < ra > 2k (m = A 
If we now use values of ra, given by the inequality (9) we find that: 


Tee Pag cient 


1 
ial == 2k (»- 3) 


In a similar manner we can show that 


Гра. == 2k (m а 1) 


$ "a2 = а= k| 
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Substituting these values of the limits in (10) we obtain a recursion formula for 
„Рои(го»). This formula must be supplemented by a formula derived in a similar fashion 


for man 24-1: 


2k(m—1) T94—1 +k 
T25—1 nci agarose) Год 1 MP 2n—2(Tan—2) 
Pe eee dra d tala заа ыз HM 
ті 2n—1 ( 2n 1) 9k (ems n 9k eee n 
lron—1—Al 2k(m—1) 


where we assume 1Р,, = Рә, in both (10) and (11). 
By repeated use of equations (10) and (11), and starting from Р„(т„) we can obtain 
solutions for any „Р, 


More convenient expressions can be obtained by repeated di fferentiation of expres- 
sions (10) and (11) with respect to. It can then be shown that: 


don Se dm Pos (lon) Ne um 2i о 

dr2n—2 ue cdd п— т] (2k) e 
d2n—3 MP on—1(Ten—1) а, atat? 2m Е 2 
ар КЕ xm пт) js п 


To evaluate a given ,P(r) using (12) or (13) we must integrate the latter expres- 
sions the appropriate number of times, the constants of integration being determined 
by applying the conditions that all the derivatives are continuous functions of r for 
values of r less than the maximum, and that all derivatives vanish at the maximum 
value of r (see Rayleigh 1919). We give below the solutions for 5 equal vectors: 


5/28 — r$ 
Руы) = ев O 
и ot 
Р) = LATTE эш Em k < rg 3k 
125k3r, — 75k2r2 + 15478 — гі 
Ру) = В E TS А oh = k 


In conclusion it should be stated that the present method, while in principle ca- 
pable of being applied to the same problems as those treated by Rayleigh and Chan- 
drasekhar, does not lead to general expressions of sufficiently simple form to allow 

application to very large numbers of vectors of known amplitude distribution. On the 
. other hand for small numbers of vectors as for example the photon--recoil problem 
mentioned in the introduction, the present method of deriving the appropriate sub- 


-ranges of r and the integrals for P(r) is much simpler to use than methods previously 
used, 


E — PRU 
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Д. Оконор, Упрощенный метод решения особого случая безпорядочного движения. 


Представлено упрощенный метод счёта вероятности нахождения в данном 
пределе амплитуды равнодействующего вектора созданного прибавлением конеч- 
ного числа векторов о известной амплитуде и случайной ориентировке. Для ма- 
лого числа векторов неодинаковой длины, этот метод проще методов известных 


_ в литературе. Дано рекуретные формулы на функции вероятности для множества 


.— векторов равной длины. 


= 
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